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ABSTRACT 


A  modified  version  of  the  Chapman-Korst  model  involving  mass  diffusion, 
equilibrium  thermochemistry,  and  appropriate  semiempirical  relations  is 
utilized  to  predict  the  overall  flow  characteristics  within  the  plume-induced 
separated  flow  region  during  synergetic  maneuvering  of  a  lifting  reentry  ve¬ 
hicle.  The  exhaust  plume  selected  in  conjunction  with  the  separated  flow 
analysis  is  predicted  on  the  basis  of  the  method-of-characteristics  solution 
with  an  average  specific-heat  ratio  of  1.38.  Thermochemical  equilibrium 
is  assumed  to  exist  within  the  plume-induced  separated  flow  region,  and 
equilibrium  temperature  curves  for  various  altitudes,  vehicle  Mach  num¬ 
bers,  and  propellant-air  mixture  ratios  are  presented.  The  strip  method 
is  incorporated  in  the  analysis  to  estimate  the  three-dimensional  separated 
flow  characteristics  caused  by  the  vehicle  geometry  and  vehicle  attitude 
during  the  maneuver.  The  extent  of  validity  of  the  strip  method  is  estab¬ 
lished  with  the  aid  of  a  cross-flow  boundary  layer  solution.  On  the  basis 
of  the  present  analytical  model,  the  propulsion  effects  on  aerodynamic 
characteristics  of  a  simple  lifting  body  configuration  with  delta  planform 
and  triangular  cross  section  under  a  wide  spectrum  of  freestream  condi¬ 
tions  have  been  evaluated  and  illustrated  by  appropriate  design  charts. 
Although  the  usefulness  of  the  prediction  method  has  been  demonstrated, 
additional  relevant  data  are  required  to  finalize  the  theory  development 
and  to  ensure  a  broad  applicability  of  the  analytical  model. 

(This  abstract  is  subject  to  special  export  controls  and  each  transmittal 
to  foreign  governments  or  foreign  nationals  may  be  made  only  with  prior 
approval  of  the  Air  Force  Flight  Dynamics  Laboratory.) 
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Section  1 
INTRODUCTION 

When  a  lifting  reentry  vehicle  (LRV)  undertakes  aerodynamic  maneuvering  to  effect 
orbital  plane  changes  (synergetic  maneuvering),  two  modes  of  operation  are  possible: 
aeroglide  turn  and  aerocruise  turn.  Aeroglide  turn  refers  to  the  maneuver  using  the 
aerodynamic  turning  capability  of  the  vehicle  in  a  normal  unpowered  glide;  propulsion 
is  required  only  for  the  initial  deorbit  impulse  and  for  achieving  the  desired  orbit  or 
flight  condition  following  completion  of  the  turning  maneuver.  Aerocruise  turn  refers 
to  the  orbital  change  maneuver  achieved  by  varying  the  vehicle  aerodynamic  forces 
through  propulsive  thrust  manipulation.  From  the  viewpoint  of  heating  restraints,  the 
aerocruise  maneuver  is  superior  to  its  aeroglide  counterpart  for  synergetic  plane 
changes  (Ref.  1). 

However,  the  presence  of  propulsive  thrust,  and  an  associated  jet  plume  issuing  from 
the  rocket  nozzle  during  the  aerocruise  maneuver  can  produce  surface  flow  separation 
and,  consequently,  undesirable  effects  on  vehicle  aerodynamics.  This  separated  flow 
phenomenon  occurs  because  at  high  altitudes  the  expanded  rocket  exhaust  plume  would 
resemble  a  large  hemisphere -cylinder  afterbody  which  may  induce  an  adverse  pressure 
gradient  too  great  for  the  boundary  layer  to  remain  attached  on  the  vehicle  surface. 
Experimental  studies  conducted  by  NASA  (Refs.  2  through  5)  and  Mithras  (Refs.  G  and  7) 
have  revealed  that  jet  pluming  can  indeed  incur  flow  separation  on  the  body  surface 
under  simulated  high  altitude  flight  conditions. 

For  sufficiently  extensive  plume-induced  flow  separation,  the  attendant  local  pressure 
elevation  could  result  in  significant  losses  in  lift,  stability,  and  control  of  the  lifting 
reentry  vehicle.  Such  losses  in  performance  would  seriously  reduce  the  synergetic 
maneuvering  potential  of  the  vehicle,  and  the  high  L/D  features  brought  forth  by  the 
aerocruise  maneuver  would  be  greatly  compromised. 

An  efficient  design  of  high-performance  lifting  reeentry  vehicles  thus  entails  configura¬ 
tion  development  as  well  as  a  careful  study  of  the  exhaust  plume  effects  on  vehicle 
aerodynamic  characteristics  under  realistic  operating  conditions.  The  present  study 
is  concerned  exclusively  with  the  latter  aspect  and  encompasses  the  following  objectives: 

•  Development  of  an  analytical  scheme  to  predict  the  overall  flow  characteristics 
in  the  plume -induced  separation  region  on  lifting  reentry  vehicles  for  a  wide 
spectrum  of  flight  variables 

•  Estimation  of  the  effects  of  plume-induced  separation  on  the  aerodynamic 
characteristics  of  lifting  surfaces 

•  Identification  of  important  physical  parameters  which  require  further  experi¬ 
mental  information  so  that  a  meaningful  and  proficient  wind  tunnel  test  pro¬ 
gram  for  the  follow-on  study  can  be  formulated 
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Specific  tasks  and  various  disciplines  involved  in  the  theory  development  are  delineated 
in  the  text.  For  convenience  and  for  order  of  discussion,  these  tasks  are  classified  in 
terms  of  several  flow  regions,  each  essentially  embracing  a  particular  discipline. 
Figure  1  depicts  these  flow  regions  along  with  the  analytical  scheme,  in  flow -diagram 
form,  for  predicting  the  overall  propulsion  effects  on  LRV  aerodynamic  characteristics. 
Basic  assumptions  and  salient  features  associated  with  each  flow  region  are  listed  in 
Table  I  to  provide  a  ciear  perspective  of  the  analytical  tasks  performed  during  this  study. 

For  completeness,  detailed  derivations  of  the  pseudo-two-dimensional  separated  flow 
theory  and  the  cross-flow  perturbation  theory  developed  during  the  study  are  included 
in  the  Appendixes.  Also  reviewed  are  the  degree  of  approximation  and  limitations  of 
the  prediction  method  as  a  result  of  the  lack  of  relevant  data.  Finally,  the  analytical 
model  is  demonstrated  by  means  of  a  detailed  analysis  of  a  simple  lifting  body  configu¬ 
ration  witn  a  delta  planform  and  triangular  cross  section  for  a  wide  range  of  flight 
conditions.  Appropriate  design  charts  and  easy-to-use  semiempirical  formulas  are 
included. 
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Table  I  Summary  of  Problem  So 
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e  Exponential-polynomial  formula  for  Chapman 
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ANALYSIS  AND  DISCUSSION 


This  section  presents  a  description  of  the  analytical  model  for  plume-induced  flow 
separation.  This  model  asserts  that  when  the  flow  field  on  a  LRV  surface  encounters 
a  sufficiently  strong  adverse  pressure  gradient  induced  by  the  exhaust  plume,  a 
constant-pressure  shear  layer  will  emanate  from  some  point  on  the  vehicle  surface 
and  intersect  with  the  corresponding  shear  layer  along  the  plume  boundary.  A  sepa¬ 
rated  flow  region  and  a  trailing  shock  system  are  thus  formed,  accompanied  by  varia¬ 
tion  of  local  flow  properties.  Important  flow  regions,  related  phenomena,  and  the 
main  tasks  involved  are  listed  in  the  Introduction.  In  this  section,  the  discussion  is 
devoted  to  detailed  descriptions  of  these  different  flow  regions,  with  emphasis  on 
existing  theories  and  unique  features  associated  with  this  class  of  flow  separation  as 
well  as  simplifications  and  refinements  underlying  the  present  analytical  approach. 

1.  CANDIDATE  PROPULSION  SYSTEMS 

Potential  rocket  engine  candidates  for  representative  lifting  reentry  vehicles  were 
reviewed  to  determine  typical  propulsion  system  parameters.  While  the  present  study 
is  not  directed  specifically  to  the  analysis  of  a  particular  vehicle  design,  this  approach 
was  selected  in  order  to  limit  the  range  of  values  to  those  which  would  bracket  candidate 
systems.  From  this  review,  two  propulsion  systems  have  emerged  as  prime  candidates , 
owing  to  their  superior  engine  performance:  hydrogen -fluorine  (Ho”Fo)  ^  hydrogen- 
oxygen  (H„-0„).  The  important  propulsion  parameters  for  the  hydrogen-fluorine  sys¬ 
tem  analyzed  during  a  concurrent  USAF  study  (Ref  8)  and  a  similar  hydrogen-oxygen 
system  are  presented  in  Tables  n  and  III.  Columns  1  and  2  in  these  tables  present 
values  for  a  small  "throttleable"  bell  nozzle;  and  columns  3,  4,  and  5  present  values 
for  a  larger  bell  nozzle  engine  equivalent  in  rated  thrust  to  the  small  engine.  The  two 
engines  (large  and  small  bells)  can  produce  identical  thrust  for  synergetic  maneuvers 
as  shown  in  Columns  2  and  3,  although  for  considerably  different  chamber  pressures 
and  exit  diameters.  Variation  of  oxidizer-to-fuel  ratio  (O/F)  for  the  particular  oxidizer 
and  fuel  combinations  aie  not  expected  to  be  significant  within  the  limits  of  the  present 
study.  The  remaining  terms  are  the  primary  parameters  which  are  to  be  varied  over  the 
indicated  excursion  range. 

In  the  present  study  which  is  concerned  with  exhaust  plume  effects,  a  single  set  of  repre¬ 
sentative  propulsion  parameters  for  both  systems  was  selected  for  detailed  analysis. 

The  parameters  include:  chamber  pressure  of  500  psi,  area  ratio  of  60,  and  nozzle  half¬ 
angle  of  9. 5  degrees  with  O/F  =  12  and  6  for  the  hydrogen-fluorine  system  and  the 
hydrogen -oxygen  system,  respectively.  Emphasis  on  exhaust  plume  description  also 
attaches  prominence  to  the  prediction  of  propellant  thermochemical  properties  discussed 
in  the  next  section. 


Table  n  Propulsion  Parameters  (H9-F9  System) 


2.  THERMOCHEMISTRY  IN  NOZZLE  AND  IN  PLUME  FLOW  FIELD 


a.  Propellant  Combustion  and  Nozzle  Flow 

The  overall  thermodynamic  properties  in  the  nozzle  and  the  exhaust  plume  flow  fields 
associated  with  the  candidate  propulsion  systems  can  be  predicted  by  coupling  thermo¬ 
chemistry  with  fluid  flow  analysis.  In  the  combustion  chamber,  where  the  fluid  motion 
is  slow,  equilibrium  chemical  calculations  can  be  performed  to  ascertain  the  basic 
thermochemical  parameters,  such  as  gas  temperature,  species  concentration,  and 
specific  heat  ratio,  Y  .  In  recent  years,  calculation  schemes  to  determine  these  quan¬ 
tities  have  been  standardized  and  continuously  improved  by  NASA-Lewis.  As  a  result, 
a  versatile  equilibrium  thermochemical  computer  program  is  now  available  for  universal 
use  (Ref  9).  A  complete  listing  of  the  NASA-Lewis  program  and  the  proper  input  format 
are  included  in  Reference  9,  which  can  be  duplicated  by  anyone  interested  in  full  use  of 
this  program.  In  this  report,  however,  calculated  thermochemical  results  for  the  can¬ 
didate  propellant  systems  are  presented  in  design  chart  form  to  permit  temperature 
interpolation  in  connection  with  separated  flow  computations. 

During  nozzle  expansion,  nonequilibrium  reactions  prevail  over  a  portion  of  the  nozzle 
for  which  tedious  computations  and  uncertain  information  on  important  reaction  rates 
are  frequently  encountered.  An  expedient  alternative  is  to  implement  Bray's  "sudden 
freezing"  criterion  (Ref  10),  in  that  the  final  flow  is  assumed  to  be  initially  character¬ 
ized  by  "shifting"  equilibrium  (infinite  reaction  rates),  but  at  some  distance  downstream 
of  the  throat,  a  "composition  freeze"  takes  place  suddenly.  This  model  essentially 
postulates  an  infinitely  thin  nonequilibrium  or  transition  zone  and  is  valid  if  three-body 
recombination  constitutes  the  controlling  reaction  process.  By  circumventing  the  non¬ 
equilibrium  consideration,  the  relatively  easy  equilibrium  or  frozen  flow  calculations 
can  be  carried  out  to  determine  the  thermochemical  quantities  in  the  chamber  and  in 
the  nozzle.  In  most  applications,  sudden  freezing  leads  to  reasonable  thermochemical 
estimates  and  occurs  at  area  ratio  less  than  5.  For  the  candidate  propulsion  systems 
considered  in  this  analysis,  the  large  area  ratio  (40  to  150)  implies  that  species  con¬ 
centrations  will  have  been  frozen  long  before  the  flow  reaches  the  nozzle  exit. 

In  this  study,  the  thermochemical  problem  in  connection  with  the  nozzle  flow  for  a 
lifting  reentry  vehicle  is  treated  as  follows: 

•  In  the  combustion  chamber,  equilibrium  chemical  computations  based  on  the 
NASA-Lewis  program  are  performed  to  determine  the  thermochemical 
properties. 

•  In  the  nozzle  portion,  a  one-dimensional  inviscid  flow  field  is  assumed  to 
exist.  For  most  applications,  low  critical  area  ratios  for  sudden  freezing 
are  observed;  hence,  the  assumption  of  composition  freezing  in  the  chamber 
is  expected  to  result  in  valid  first  approximation.  For  simplicity,  this 
chamber  freezing  model  is  postulated  along  with  a  constant  Y  value  of  1.38 
as  an  average  "frozen"  value.  Selection  of  this  relatively  large  Y  value 

is  consistent  with  the  continual  increase  in  y  and  decrease  in  c  during 
nozzle  expansion  under  frozen  flow  conditions.  The  wall  boundary  layer 
effect  is  not  expected  to  affect  the  nozzle  flow  significantly;  its  effect  on 
plume  properties  will  be  discussed  later. 


Table  IV  presents  typical  thermochemical  results  for  the  candidate  hydrogen-fluorine 
and  hydrogen -oxygen  systems,  including  equilibrium  and  frozen  temperatures  at  dif¬ 
ferent  pressure  ratios  (or  area  ratios)  that  correspond  to  different  nozzle  locations. 

Once  the  nozzle  exit  thermochemical  properties  are  properly  defined,  the  entire  ex¬ 
haust  plume  flow  field  can  be  described  by  means  of  inviscid  flow  analysis  and  gas- 
dynamic-thermochemical  coupling. 

b.  Thermochemistry  in  Plume  Flow  Field 

Beyond  the  nozzle  exit,  the  thermochemical  effects  are  of  importance  to  LRV  applica¬ 
tions  in  two  aspects: 

•  Plume  geometry  and  other  physical  features  are  influenced  by  the  thermo¬ 
chemical  properties  of  the  flow  field,  e.  g. ,  the  specific  heat  ratio,  y  , 
thereby  affecting  the  plume-induced  separation  mechanism. 

•  Diffusion  of  the  generally  fuel-rich  exhaust  products  may  occur  in  the  shear 
layer  as  well  as  in  the  base  and  plume-induced  separation  region.  The 
ensuing  mixing  between  the  exhaust  gas  and  the  ambient  air  may  initiate 
chemical  reactions,  causing  an  elevated  base  temperature*  and  possible 
overheating  of  the  vehicle  surface.  Deleterious  by-products  arising  from 
chemical  reactions,  if  existing,  may  lead  to  surface  corrosion  after  pro¬ 
longed  exposure. 

In  view  of  the  constant  7  model  for  plume  description  adopted  in  the  present  analysis, 
this  section  focuses  attention  on  establishing  conditions  where  afterburning  occurs 
and  preparing  appropriate  temperature  charts  useful  to  plume-induced  separated  flow 
predictions. 

As  already  mentioned  in  the  Introduction  and  as  will  be  further  elaborated  in  a  later 
section,  die  generalized  Crocco  relation  is  rendered  valid  as  a  result  of  the  assump¬ 
tions  of  unit  Prandtl  and  Schmidt  numbers  and  equilibrium  chemical  reactions.  The 
local  temperature  and  local  exhaust -gas/ambient-air  mixture  ratio  within  the  shear 
layer,  including  the  plume-induced  separation  region,  can  then  be  ascertained. 

Based  on  the  "flame  sheet"  concept,  the  afterburning  properties  along  any  shear  layer 
streamlines  can  be  estimated  also.  However,  from  the  results  of  Reference  11,  the 
reaction  time  and  the  ignition  delay  time  vary  inversely  with  ambient  pressure,  which 
render  afterburning  in  a  high-speed  stream  or  shear  layer  almost  impossible  at  alti¬ 
tudes  near  100,000  feet.  On  the  other  hand,  the  plume-induced  separation  problem 
is  an  exception.  In  this  case,  the  base  region  is  characterized  by  a  "dead-air"  zone 
where  the  exhaust  gas  and  the  ambient  air  will  diffuse,  mix,  and  become  stagnated, 
therefore,  a  short  characteristic  ignition  length  is  ensured  despite  the  long  reaction 
and  ignition  times,  and  the  chemical  equilibrium  model  is  generally  valid. 


♦Base  temperature  and  pressure  refer  to  either  base  or  plume-induced  separation 
region. 


Table  IV  Equilibrium  Thermochemical  Calculations  for  H  -F~  and  H0-0_  Systems 
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For  the  present  LRV  application,  only  the  dead-air  case  is  considered  for  afterburning 
calculations.  The  equilibrium  model  is  adopted  and  design  curves  are  prepared  which 
correlate  equilibrium  base  temperature  with  exhaust-gas  mixture  ratio,  CR  (ratio  of 
exhaust  mass  to  total  mass  of  exhaust  gas  and  air),  for  various  ambient  pressures 
(or  altitudes)  and  equivalent  heat  of  formation,  A  hx.  The  latter  quantity  is  defined  as 
the  enthalpy  in  calories  per  mole  of  air  converted  from  the  kinetic  energy  of  the  high¬ 
speed  stream.  In  other  words,  when  the  free  stream  is  stagnated  in  a  dead-air  region, 
the  chemical  system  within  this  region  will  experience  an  increase  in  internal  energy 
arising  from  the  kinetic  energy  conversion.  Thermochemical  calculations  are  there¬ 
fore  performed  by  incorporating  this  heat  of  formation  term  which  depends  on  Mach 
number  and  altitudes,  as  shown  in  Figure  2.  Figures  3  through  19  depict  the  appro¬ 
priate  temperature  design  curves.  It  should  be  remarked  that  the  temperature  and 
CR  values  on  these  graphs  are  synonymous  to  T~  and  Y~  used  in  conjunction  with 
the  plume-induced  separated  flow  iteration  scheme  described  in  Section  5  and  Appendix  I. 

3.  EXHAUST  PLUME  DESCRIPTION 

For  a  balanced  exhaust  jet  (nozzle  exit  pressure  equal  to  ambient  pressure),  the  flow 
field  can  be  successfully  described  by  means  of  a  viscous  shear  layer  model  shown  in 
Figure  20a.  Analytical  solutions  have  been  obtained  not  only  for  laminar  jets  but  also 
for  turbulent  jets  with  diffusion  and  chemical  reactions  (Refs  12  to  15).  However,  in 
most  high-altitude  applications,  the  rocket  exhaust  plume  is  under-expanded  (nozzle 
exit  pressure  greater  than  ambient  pressure)  for  which  its  basic  structure  is  illustrated 
in  Figure  20b.  A  comparatively  simple  method  of  treating  undex expanded  plume  prob¬ 
lems  is  to  combine  the  solution  of  the  invlscid  exhaust  plume  with  that  of  a  viscous  mix¬ 
ing  layer  (shear  layer)  which  is  assumed  to  grow  continuously  along  the  plume  boundary. 
The  near-field  region,  which  features  an  extensive  inviscid  flow  field  and  a  thin  shear 
layer,  is  particularly  amenable  to  this  approach.  Inasmuch  as  plume-induced  flow 
separation  extends  within  the  near-field  region  only,  the  phenomena  of  Mach  disc,  shock 
reflection,  and  turbulent  wake  occurring  in  the  intermediate  and  far-field  regions  are 
neglected  in  the  present  analysis.  Only  the  inviscid  portion  of  the  exhaust  plume  is  con¬ 
sidered  here;  the  shear  layer  problem  is  discussed  subsequently. 

a.  Inviscid  Plume  Prediction 

The  inviscid  plume  contour  and  its  Internal  properties  can  be  described  accurately  by 
employing  the  method  of  characteristics.  The  basic  principles  underlying  the  method- 
of-characteristics  (MOC)  solution  for  axisymmetric  exhaust  plumes  and  formulation  of 
a  computation  scheme  are  discussed  in  detail  in  Reference  16.  To  oroceed  with  the 
computation,  the  one -dimensional  or  axisymmetric  MOC  solution  for  the  nozzle  flow 
is  first  obtained  to  establish  the  nozzle  exit  conditions  which  represent  the  initial  con¬ 
ditions  for  the  exhaust  plume.  Then,  step-by-step  MOC  numerical  computation  is 
carried  out  for  the  plume  flow  field  until  the  free  boundary  or  ambient  pressure  con¬ 
dition  is  satisfied.  One  computation  scheme  is  to  construct  the  left-running  charac¬ 
teristics  from  point  to  point  over  the  entire  flow  field  to  determine  the  local  flow 
properties,  including  the  locus  for  the  ambient  pressure  points  or  the  free  boundary. 

A  boundary  shock  is  also  formed  downstream  of  the  nozzle  exit  as  a  result  of  coalesc¬ 
ing  compression  waves  reflected  from  the  boundary. 
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Figure  4  Hydrogen-Fluorine/Air  Equilibrium  Temperatures  for  Ah^  =  5000  Cal/Mole 
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Figure  7  Hydrogen- Fluorine/Air  Equilibrium  Temperatures  for  Ah^  =50,000  Cal/Mole 
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Figure  9  Hydrogen- Fluorine /Air  Equilibrium  Temperatures  forAhf  =  100,000  Cal/Mole 


Figure  11  Hydrogen-Oxygen/Air  Equilibrium  Temperatures  for  Ah^  =  0 
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Figure  12  Hydrogen-Oxygen/Air  Equilibrium  Temperatures  for  =  86,000  Cal/Mole 
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Figure  13  Hydrogen-Oxygen/ Air  Equilibrium  Temperatures  for 
hf  =  149,000  Cal/Mole 
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Figure  15  Liquid  Hydrogen/ Liquid  Fluorine/Air  Equilibrium  Temperatures  for 

O/F  =  12,  P  =  500  psi 
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Figure  17  Liquid  Hydrogen/ Liquid  Oxygen/Air  Equilibrium  Temperatures  for 
O/F  =6,  Pc  =  500  psi 
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Figure  18  Liquid  Hydrogen/Liquid  Oxygen/ Air  Equilibrium  Temperatures  for 

O/F  =  6,  P  =  500  psi 
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Figure  20  Sketches  of  Plume  Structure 
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In  general,  two  types  of  MOC  problems  in  terms  of  boundary  conditions  can  be 
classified: 


•  Quiescent  ambient  condition.  The  only  boundary  condition  to  be  satisfied  is 
that  the  ambient  pressure  must  be  matched  at  the  free  boundary.  Furthermore, 
according  to  the  MOC  model,  all  other  thermodynamic  quantities  are  constant 
along  the  plume  contour.  Figure  21  illustrates  a  typical  plume  boundary  under 
quiescent  ambient  condition  and  several  constant-Mach  contours  in  the  interior. 
Figure  22  depicts  several  important  plume  boundary  parameters  versus  ambient 
pressure,  each  being  constant  along  the  entire  plume  contour  for  a  given  ambient 
pressure. 

•  High-speed  freestream  condition.  Because  of  external  stream  impinging  on  the 
plume,  the  thermodynamic  properties  along  the  plume  boundary  will  vary  from 
point  to  point.  The  "local"  plume  contour  can  be  determined  iteratively  with 
the  aid  of  the  oblique  shock  or  Newtonian  concept  with  which  the  boundary 
pressure  of  the  characteristic  mesh  must  match  the  impact  pressure.  Figure  23 
shows  typical  contours  under  high-speed  freestream  conditions. 

A  computer  program  based  on  this  MOC  solution  currently  in  use  at  LMSC  has  been 
proven  efficient  and  commands  a  wide  applicability.  Satisfactory  comparisons  with  air 
plume  data  for  a  number  of  sample  problems  have  also  been  demonstrated  (Refs  17  and 
18).  In  addition  to  frozen  flow  and  constant  Y  ,  the  program  also  accommodates  equili¬ 
brium  chemistry  in  each  characteristic  mesh  as  well  as  an  option  to  estimate  the  bound¬ 
ary  layer  effects  during  nozzle  expansion.  However,  chemical  kinetics,  phase  transition, 
and  three-dimensional  plume  effects  cannot  be  accurately  predicted  with  the  LMSC  MOC 
plume  program  in  its  present  version. 

While  the  MOC  solution  has  been  extensively  applied  during  this  study,  relevant  design 
charts  and  simplified  formulas  for  exhaust  plume  description  have  been  developed  so 
that  important  plume  properties  for  the  candidate  propellant  systems  can  be  determined 
quickly.  In  this  respect,  the  following  two  postulates  were  made: 

•  Consistent  with  the  dead-air  hypothesis,  the  plume  surface  enshrouded  by  the 
separated  flow  region  can  be  regarded  as  being  exposed  to  zero  ambient  air 
velocity.  As  this  plume  portion  is  pertinent  to  the  trailing  shock  determination 
and  the  separated  flow  solution  (to  be  discussed  in  a  later  section),  design  curves 
are  constructed  only  on  the  basis  of  the  quiescent  environment. 

•  From  the  discussion  on  thermochemistry  in  the  preceding  section,  it  is  quite 
certain  that  forzen  chemical  composition  would  be  established  within  a  short 
distance  from  the  nozzle  throat.  However,  because  of  continuous  flow  expan¬ 
sion,  the  V  value  would  increase  with  decreasing  temperature.  A  large  y 
value  in  the  plume  flow  field  is  thus  anticipated  (likely  to  be  of  the  order  1.4) 
even  if  the  equilibrium  y  value  were  as  low  as  1. 1  in  the  combustion  chamber. 
As  already  indicated  in  the  preceding  section,  a  constant  y  of  1.38  is  selected 
in  the  present  analysis.  The  typical  effect  of  y  on  plume  contour  is  depicted 
in  Figure  24. 


CONSTANT  PROPERTY  CONTOUR 


Figure  21  MOC  Jet  Plume  in  Quiescent  Air 
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Figure  23  MOC  Jet  Plume  at 


With  the  plume  model  selected,  the  exhaust  plume  properties  for  hydrogen-fluorine 
and  hydrogen -oxygen  mixtures  are  obtained  by  carrying  out  the  MOC  computations 
together  with  the  following  assumptions: 

•  Axisymmetric  with  zero  angle  of  attack 

•  Area  ratio  =  60 

•  Ambient  air  velocity  =  ,0 

•  Chamber  pressure  =  500  psi 

•  Chamber  temperature  =  8050°R 

•  V  =  1.38 

•  Nozzle  half-angle  =  9.5  degrees 

•  Frozen  species  concentrations 

By  selecting  a  particular  V  ,  the  propulsion  system  is  not  explicitly  identifiable  except 
by  implication  of  the  chamber  temperature  value.  However,  the  exhaust  plume  contour 
is  affected  by  the  chamber  temperature  primarily  through  the  nozzle  exit  Mach  number. 
According  to  calculated  results  presented  in  Ref  19,  this  effect  is  not  a  strong  one. 
Furthermore,  the  assumed  chamber  pressure  of  500  psi  would  lead  to  no  loss  of  gen¬ 
erality  since  the  design  curves  were  constructed  in  terms  of  jet-to-ambient  pressure 
ratio,  as  depicted  in  Figure  25. 

For  development  of  simplified  prediction  formulas,  two  methods  are  proposed  herein: 
(1)  Formulas  to  predict  five  boundary  points  based  on  the  MOC  solution,  and  (2)  The 
Latvala -Anderson  model  to  predict  plume  contour  in  the  vicinity  of  the  nozzle  exit. 

(1)  The  Five- Point  Method.  On  the  basis  of  the  MOC  solutions,  prediction 
formulas  have  also  been  developed  which  permit  estimates  of  plume 
contours  involving  the  effects  of  freestream  Mach  number,  area  ratio, 
nozzle  half-angle,  etc. 
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Figure  25  Exhaust  Plume  Contours  for  Various  Ambient  Pressures 
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The  corresponding  prediction  curves  are  presented  in  Figures  26 
through  28.  It  should  be  mentioned  that  these  formulas  are  valid  for 
non-zero  freestream  Mach  number  also,  but  they  are  based  on  y  =  1.38 
and  zero  angle  of  attack  for  the  plume. 


(2)  The  Latvala -Anderson  Method.  The  Latvala-Anderson  Method  (Refs  20 
and  21)  provides  a  simplified  scheme  to  predict  the  initial  portion  of  a 
highly  underexpanded  jet  boundary.  The  basic  hypothesis  underlying 
the  method  is  that  the  jet  boundary  is  composed  of  a  series  of  connecting 
circular  arcs,  each  being  describable  by  the  following  "length"  formulas 
(Fig.  29) 


rn+ 1  _ 

■  (1  +  cos  an  +  I)  (£)n+1' 

1/2 

rex 

-  <*  +  cos  V  (jF-)n  - 

*x/rex  . 

/  rn+l^rex  A  t  / 

AQ 

rn^rex 

I  v^r  7  ( “ 

n  "  2 

(8) 


(9) 


The  initial  wave  angle  is  determined  from  the  Prandtl- Meyer  expression 


where  v  is  the  Prandtl-Meyer  angle,  and  M  is  the  average  Mach  number. 
The  corresponding  turning  angle  is  then  determined  from 
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Figure  28  Local  Plume  Radius  as  a  Function  of  Maximum  Plume  Radius 
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Continuing  the  forward  integration  using  local  Prandtl-Meyer  angle,  average 
Mach  number,  and  isentropic  area  ratio,  the  entire  spherical-arc  jet  bound¬ 
ary  surface  for  the  fore  part  of  the  plume  can  be  ascertained.  Figure  30 
depicts  the  comparison  between  the  Latvala -Anderson  solution  and  the  MOC 
solution  for  several  typical  air  plumes  which  show  excellent  agreement  for 
the  region  investigated. 

For  further  approximation,  graphical  construction  of  the  circular-arc  initial 
jet  piume  boundary  can  be  made.  Figure  31  shows  the  variation  of  the  radius 
as  a  function  of  nozzle  Mach  number,  Mex,  for  an  air  plume.  The  curve  can 
also  be  represented  by  the  following  simple  formula: 


_  15.  25  +  i  70 

r  M  -  0  377  13  78 

ex  ex 


(12) 


Although  the  Latvala- Anderson  method  is  merited  with  simplicity  in  execution, 
its  one-dimensional  model  admits  only  an  average  Mach  number  and  other 
average  flow  quantities  across  the  plume  for  which  the  pressure  boundary 
condition  is  not  satisfied.  Hence,  aside  from  the  plume  contour  near  the 
nozzle  exit,  the  important  plume  parameters  along  the  jet  boundary  cannot 
be  determined  by  this  method  Additional  information,  such  as  the  design 
curves  presented  in  Figure  22,  must  be  supplemented  in  order  to  pursue  the 
plume-induced  separated  flow  prediction. 

b.  y  Variation  in  Exhaust  Plume 

As  mentioned  previously,  rapid  high-altitude  plume  expansion  under  frozen  flow  condi¬ 
tions  would  result  in  a  decrease  in  gas  temperature  and  a  corresponding  rise  in  the  Y 
value.  A  rigorous  approach  to  plume  prediction  entails  a  continuous  adjustment  of  the 
Yvalue  in  the  characteristic  mesh.  On  the  other  hand,  a  simplified  model,  as  adopted 
herein,  is  to  select  a  constant  Y corresponding  to  anticipated  boundary  temperatures 
Although  somewhat  arbitrary,  a  constant  Y  of  1. 38  has  been  chosen  on  the  basis  of  ob¬ 
serving  a  number  of  MOC  plumes  with  various  Y  values  for  the  nozzle  exit  conditions. 
Evidently,  the  Latvala-Anderson  method  with  Y  =  1.4  would  lead  to  practically  the 
same  answer. 

c.  Nozzle  Boundary  Layer  Effects 

The  growth  of  boundary  layer  in  the  rocket  nozzle  affects  the  exhaust  plume  properties 
in  two  ways: 

•  Modification  of  the  nozzle  shape  and  area  ratio  due  to  the  growth  of  displace¬ 
ment  thickness  along  the  wall 

•  Change  in  plume  contour  due  to  the  presence  of  low-speed  boundary  layer  at 
the  nozzle  exit 
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(Reproduced  from  Reference  20) 
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Figure  30  Jet  Boundaries  Calculated  by  Approximate  Method  Compared 
With  Method  of  Characteristics  Solution  from  NASA  TRR-6 
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The  first  problem  has  been  investigated  quite  extensively,  particularly  by  Bartz  et  al 
(Refs.  22  and  23)  in  conjunction  with  the  evaluation  of  rocket  nozzle  heat  transfer 
characteristics  and  possible  degradation  of  performance  as  a  result  of  boundary  layer 
formation.  By  employing  the  turbulent  nozzle  flow  model  discussed  in  Reference  22, 
the  displacement  thickness  at  the  nozzle  exit  for  the  hydrogen-fluorine  plume  considered 
here  (AR=60)  was  found  to  be  approximately  2  percent  of  the  exit  radius.  With  this 
magnitude  of  displacement  thickness  growth,  the  nozzle  geometry  will  be  modified  only 
slightly,  for  which  any  significant  change  in  plume  configuration  is  not  expected  to  be 
incurred. 

However,  upon  leaving  the  nozzle  exit,  the  low-speed  wall  boundary  layer  must  undergo 
greater  expansion  than  its  high-speed  counterpart  in  the  center  core  so  as  to  satisfy 
the  pressure  boundary  condition.  As  a  result,  the  plume  boundary  would  become  en¬ 
larged,  particularly  in  the  near-field  region  (Ref.  24).  To  assess  this  effect,  the 
boundary  layer  option  of  the  LMSC  MOC  plume  program  was  utilized  which  calculates 
the  boundary  layer  profile  in  the  nozzle  and,  using  M=1  as  the  lowest  Mach  number, 
subsequently  computes  the  Prandtl-Meyer  expansion  at  the  nozzle  exit  corner.  The 
program  also  computes  the  exit  Reynolds  number  and  selects  laminar  or  turbulent 
profile  accordingly. 

Using  the  hydrogen-fluorine  system  once  again,  the  plume  contours  showing  the  bound¬ 
ary  layer  effects  were  obtained  (Fig.  32).  While  a  laminar  profile  has  been  predicted 
based  on  a  transitional  Reynolds  number  of  106,  the  turbulent  case  is  also  included, 
since  in  rocket  nozzles,  the  flow  is  almost  always  turbulent.  The  comparison  presented 
in  Figure  32  shows  a  considerable  enlargement  in  plume  size  when  the  boundary  layer 
effect  is  considered.  Other  parameters,  such  as  Internal  shock  location,  boundary 
Mach  number,  etc. ,  will  also  change  accordingly.  These  changes  will  strongly  affect 
the  ensuing  plume-induced  separation  prediction  and  must  be  understood  thoroughly. 
Discrepancy  between  inviscid  plume  prediction  and  schlieren  observation  for  selected 
cases  obtained  from  recent  NASA-Langley  plume-induced  separation  study  (Ref.  25) 
seems  to  indicate  the  boundary  layer  effect  qualitatively.  However,  data  available  are 
not  sufficiently  copious  to  draw  a  final  conclusion  at  this  time. 

d.  Effect  of  Sting  Support  Concentric  to  the  Nozzle 

Some  attention  was  accorded  to  this  problem  because  a  sting  support  concentric  to  the 
jet  nozzle  was  used  in  all  the  NASA-Langley  tests  (Ref.  25).  By  considering  the  same 
area  ratio  and  nozzle  angle,  it  was  found  that  the  sting  would  not  yield  any  effect  on 
the  inviscid  plume  properties  (Fig.  32). 

e.  Three-Dimensional  Plume  Effect 

Although  only  an  axisymmetric  nozzle  is  considered  in  this  analysis,  the  exhaust  plume 
will  become  three-dimensional  when  the  vehicle  is  at  angle  of  attack  or  yaw.  A  rigorous 
analysis  for  describing  three-dimensional  plumes  requires  tedious  numerical  computa¬ 
tions.  However,  for  a  reasonable  estimate  of  the  three-dimensional  plume  character¬ 
istics  which  are  coupled  with  the  separated  flow  analysis,  a  simple  patching  of  solutions 
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would  suffice.  For  example,  the  flow  fields  on  the  windward  side  and  the  leeward  side 
of  a  pitching  lifting  reentry  vehicle  can  be  treated  separately,  each  being  exposed  to 
its  own  local  surrounding.  The  exhaust  plume  properties  in  each  main  region  are, 
therefore,  governed  by  the  local  freestream  conditions.  On  the  compression  side,  the 
plume  would  become  more  slender,  whereas  on  the  expansion  side  the  plume  would  be 
enlarged.  The  validity  of  this  model  encompassing  these  two  distinct  regions  is  shown 
in  Figure  33  wherein  predicted  and  observed  plume  contours  at  angle  of  attack  are 
compared.  This  simple  three-dimensional  plume  model  is  chosen  in  this  analysis  in 
carrying  out  the  plume  and  separated  flow  predictions. 

4.  VISCOUS  MIXING  LAYER. 

The  inviscid-viscous  plume  model  adopted  in  the  present  analysis  is  composed  of  an 
inviscid  plume  and  a  relatively  thin  viscous  mixing  layer  (shear  layer)  superimposed 
along  the  plume  boundary  (Region  (3)  in  Fig.  1).  The  shear  layer  is  assumed  to  be 
describable  by  the  Prandtl  boundary  layer  equations  for  which  the  pressure  gradient 
across  the  shear  layer  is  zero.  In  view  of  the  hypothesis  of  an  average  pressure  within 
the  dead-air  separated  flow  region,  the  shear  layer  is  essentially  of  the  flat -plate  type 
with  zero  streamwise  pressure  gradient,  and  its  solution  is  self-similar  except  in  cases 
where  finite  boundary  layer  thickness  or  disturbance  occurs  at  the  initial  point  of  the  • 

shear  layer  (Fig.  1).  Aside  from  the  plume  boundary,  the  separation  region  induced 
by  the  enlarged  plume  is  also  characterized  by  a  shear  layer  emanating  from  the  sepa¬ 
ration  point.  In  general,  on  the  plume  side,  the  shear  layer  is  fully  developed  turbulent 
because  of  the  turbulent  boundary  layer  prevailing  in  the  nozzle.  On  the  other  hand, 
under  high-altitude  conditions  in  most  LRV  applications,  the  shear  layer  on  the  stream- 
side  is  generally  pure  laminar  (Fig.  1).  In  this  study,  the  flow  separation  and  plume 
shear  layer  are  assumed  to  be  either  pure  laminar  or  fully  developed  turbulent;  tran¬ 
sitional  shear  layer  is  excluded  entirely.  Salient  features  of  the  shear  layer  solutions 
are  discussed  in  the  following  paragraphs. 

a.  Laminar  Shear  Layer 

A  constant-pressure  laminar  shear  layer  can  be  described  by  a  set  of  partial  differential 
equations  identical  with  those  for  a  flat-plate  boundary  layer  but  with  different  boundary 
conditions.  As  discussed  by  Chapman  (Ref.  26)  the  shear-layer  momentum  equation  with 
zero  initial  boundary  layer  thickness  in  terms  of  a  similarity  velocity  function,  f,  re¬ 
duces  to  the  familiar  Blasius  equation 

t'"  +  ff"  =0  (13) 


where 
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and  the  boundary  conditions  are 
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Figure  33  Comparison  of  Predicted  and  Measured  Exhaust  Plume  Boundary  at 
Angle  of  Attack 
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f'  (oo)  =  l,  f(0)  =  0,  and  f  (-«)  =  = 


The  last  boundary  condition  indicates  that  the  average  velocity  in  the  separated  flow  or 
base  region,  ug,  is  very  low,  as  compatible  with  the  dead-air  assumption.  Numerical 
solution  to  equation  (13)  yields  the  similarity  profile  illustrated  in  Figure  34. 

To  incorporate  the  laminar  profile  in  the  separated  flow  computer  program,  it  was 
found  more  expedient  to  describe  the  profile  numerically.  Based  on  a  curve  fitting 
scheme,  the  velocity  profile  formula  is  obtained  as  follows: 


f'(Tj)  =  1  =  exp  (-0.5613  +  0.45491?  -  0.08868*)2  -  0.00106T13)  (15) 


where  V  =  -7.4  and  3.6  are  the  practical  lower  and  upper  limits  for  the  width  of  the 
shear  layer. 

b.  Turbulent  Shear  Layer 

The  shear  layer  initiated  at  the  rocket  nozzle  exit  and  growing  continuously  along  the 
plume  boundary  can  be  represented  approximately  by  a  constant-pressure  turbulent 
shear  layer.  For  simplicity,  the  Korst  model  is  adopted  to  describe  this  shear  layer 
for  which  the  governing  differential  equation  is  linearized  and  the  resulting  self-similar 
velocity  profile  can  be  represented  by  an  error  function  (Ref.  27).  If  the  shear  layer 
is  exposed  to  a  dead-air  region,  a  simple  error  function  formula  for  velocity  profile 
description  is  obtained  (Fig.  35): 


£  =4*  (1  +  erfT}  (16> 

e 


where  x,  y  are  the  intrinsic  coordinates  oriented  along  the  plume  boundary  and  a  is 
a  similarity  spreading  parameter  to  provide  a  macroscopic  description  of  the  turbulent 
shear  layer.  An  approximate  prediction  formula  for  a ,  either  for  a  compress  ible  free 
jet  (Ref.  27)  or  for  parallel  streams  (Ref.  28),  can  be  expressed  in  the  following  form: 


<J  =  12  +  2.758  M.  (17) 
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Figure  34  Velocity  Function  vs  Distance  Function  for  Constant  Pressure 
Laminar  Mixing  Layer 
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Figure  35  Error  Function  Velocity  Profile  for  Constant  Pressure 
Turbulent  Mixing  Layer 
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where  Mo  is  the  frcestream  Mach  number  for  either  a  single  stream  or,  in  the  case 
of  two  parallel  streams,  an  equivalent  single  stream.  A  summary  of  experimentally 
obtained  a  values  (Ref.  29)  is  presented  in  Figure  36,  including  a  comparison  with  the 
Maydew-Reed  semiempirical  estimates  (Ref.  30).  Data  scattering  shown  in  the  figure 
clearly  indicates  the  need  for  extensive  experimental  study  and  statistical  analysis. 

c.  Generalized  Crocco  Relation 

A  well  known  particular  solution  for  determining  the  total  enthalpy  profile  across  a 
boundary  layer  or  shear  layer  is  by  means  of  the  Crocco  relation.  The  main  hypothesis 
underlying  this  approach  is  unit  Prandtl  number  whereby  the  solution  of  the  energy 
equation  exhibits  a  linear  relation  between  the  local  total  enthalpy  and  the  local  stream- 
wise  velocity,  viz. , 


H  =  Au  +  B 


(18) 


where  A  and  B  are  constants  depending  on  the  appropriate  boundary  conditions. 
Thus,  the  relatively  complicated  energy  equation  is  now  decoupled  from  the  momentum 
equation  and  can  be  solved  by  employing  a  simple  scheme,  once  the  velocity  profile  is 
determined.  In  the  majority  of  cases,  the  unit  Prandtl  number  provides  satisfactory 
solutions,  thus  rendering  the  Crocco  relation  valid. 

By  the  same  token,  the  Crocco  relation  can  be  generalized  to  predict  the  species  con¬ 
centration  profile  across  the  boundary  or  shear  layer  if  the  corresponding  Schmidt 
number  (or  Lewis  number)  is  unity,  vix. , 


Yi  =  Au  +  b' 


(19) 


is  the  mass  concentration  for  species  i  and  A  ,  B*  are  the  integration 


In  short,  Incorporation  of  the  generalized  Crocco  relation  in  the  Chapman-Korst  model 
permits  a  complete  determination  of  the  velocity,  temperature,  and  concentration  pro¬ 
files  across  the  shear  layer  separating  the  exhaust  plume  from  the  surrounding  air. 

d.  Afterburning  in  Mixing  Layer 

By  virtue  of  the  generalized  Crocco  relation,  the  species  concentration  profile  across 
the  mixing  layer  is  defined,  which  encompasses  an  entire  range  of  mixture  ratios  of 
the  fuel-rich  exhaust  gas  and  the  ambient  air.  Combustion  or  afterburning  will  occur 
if  the  local  temperature  is  above  the  ignition  temperature  and  the  total  reaction  delay 
time  is  short.  In  the  analysis  the  flame  sheet  concept  is  invoked  for  which  local  com¬ 
bustion  is  assumed  to  take  place  in  each  layer  that  chararacterizes  a  particular  exhaust 
gas/air  miyture  ratio.  This  model  can  easily  accommodate  equilibrium  chemistry  and 
permits  construction  of  design  curves  showing  afterburning  characteristics 


where  Y. 
constants. 
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However,  in  the  shear  layer  the  energy  of  the  chemical  system  is  partially  converted 
into  kinetic  energy,  resulting  in  a  reduction  of  the  local  static  enthalpy  and  temperature. 
The  temperature  reduction  and  the  attendant  high  gas  velocity  render  the  shear  layer 
*  susceptible  to  nonequilibrium  chemical  effects.  As  a  matter  of  fact,  in  the  present 

LRV  application  involving  altitudes  above  200,000  feet,  afterburning  within  the  shear 
layer  is  not  likely  to  occur,  at  least  not  in  the  equilibrium  sense.  This  aspect  is  there¬ 
fore  not  pursued  except  for  the  base  or  plume-induced  separation  region  which  can  be 
regarded  as  a  special  mixing  layer  featuring  a  single  mixture  ratio,  zero  kinetic  energy, 
and  short  characteristic  length. 

e.  Initial  Boundary  Layer  Effects 

The  Chapman  and  Korst  profiles  presented  in  Equations  15  and  16  are  applicable  only 
to  the  case  where  zero  boundary  layer  thickness  or  flow  disturbance  prevails  at  the 
separation  point  from  which  the  shear  layer  emanates.  If  the  initial  boundary  layer  is 
finite,  then  the  self-similar  feature  of  the  Chapman-Korst  profile  will  no  longer  be 
valid.  To  emphasize  this,  the  laminar  flat-plate  case  is  considered.  In  this  case,  the 
flow  field  is  characterized  initially  by  a  self-similar  Blasius  profile  but  is  later  re¬ 
placed  by  a  different  self-similar  Chapman  profile.  The  presence  of  two  distinct 
,  velocity  profiles  violates  the  similarity  requirement  and  suggests  the  existence  of  a 

transition  region  linking  these  two  distinct  similar  flow  regions.  A  sufficient  under¬ 
standing  of  the  transition  region  thus  constitutes  the  prerequisite  for  a  reliable  assess¬ 
ment  of  the  initial  boundary  layer  effect. 

On  the  plume  side,  the  boundary  layer  thickness  at  the  nozzle  exit  (initial  point  for  the 
shear  layer)  is  considered  to  be  insignificant  for  the  present  LRV  application;  therefore, 
the  initial  boundary  layer  effect  is  excluded.  However,  on  the  stream  side,  the  local 
environment  may  be  highly  rarefied  as  a  result  of  high  altitude  and/or  large  angle  of 
attack  (leeside).  The  boundary  layer  thickness  concomitant  with  the  typical  low 
Reynolds  number  in  such  a  rarefied  gas  medium  may  be  quite  substantial;  hence,  its 
effect  must  be  considered. 

A  rigorous  approach  to  the  transition  problem  requires  a  separate  solution  (Ref.  31) 
and/or  empirical  data,  particularly  for  the  turbulent  case.  An  alternate  approach  is 
to  introduce  certain  transformation  variables  so  that  the  resulting  equivalent  shear 
layer  will  retain  the  self-similarity  features  while  Implicitly  taking  account  of  the 
initial  boundary  layer  effect.  One  such  method  is  that  based  on  the  equivalent  bleed 
concept  discussed  in  Reference  28.  With  this  model,  the  governing  integral  equations 
involving  the  initial  boundary  layer  effect  are  expressed  in  a  form  identical  with  that 
involving  mass  bleed  with  zero  initial  disturbances.  However,  as  pointed  out  in  Ref¬ 
erence  32,  a  singularity  is  present  at  the  separation  point. 

A  different  method  to  resolve  the  initial  disturbance  dilemma  is  that  based  on  the  origin 
shift  model  discussed  in  References  32  through  36.  According  to  this  model,  which  is 
adopted  in  this  analysis,  an  apparent  or  equivalent  origin  is  defined  at  a  distance 
(x  ,  y  )  upstream  of  the  actual  initial  point  satisfying  the  following  criteria  (Fig.  37); 


% 
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Figure  37  Initial  Boundary  Layer  Effect  on  Plume -Induced  Flow  Separation 


•  The  initial  boundary  layer  thickness  at  the  apparent  origin  is  zero 

•  The  mass  and  momentum  fluxes  are  the  same  as  those  at  the  actual  initial 

►  point. 

Predicted  velocity  profiles  based  on  the  origin  shift  method  were  found  to  agree  quite 
satisfactorily  with  data,  whereas  those  based  on  the  equivalent  bleed  model  have  been 
less  successful  (Ref.  35). 

The  origin-shift  model  is  applicable  to  laminar  and  turbulent  shear  layers  alike,  al¬ 
though  it  was  developed  originally  for  turbulent  flows.  Recently,  this  method  has  also 
been  applied  to  the  two-stream  mixing  problem  (Ref.  35). 

Applying  the  total  mass  and  momentum  flux  balances  between  the  actual  and  the  apparent 
origins,  the  origin-shift  coordinates  are  found  to  be  (Fig.  37): 

laminar  case 
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where  all  these  symbols  are  defined  in  the  Nomenclature  section.  Since  the  initial 
disturbance  problem  is  encountered  primarily  on  the  stream  side  in  LRV  applications, 
only  the  laminar  Equations  20  and  21  are  used  in  this  analysis.  After  some  manipula¬ 
tion,  Equation  20  can  be  simplified  as  follows: 


(24) 
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where  the  constant  value  26  is  based  on  n  -  3.6  for  the  shear  layer  thickness  above 
the  plume  boundary  (Fig.  34). 

Equations  (21)  and  (24)  represent  the  final  form  used  in  conjunction  with  the  separated 
flow  analysis  in  this  report.  The  local  boundary  layer  thickness  is  determined  by 
means  of  the  approximate  formula  presented  in  Reference  37 
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where  subscript  a  represents  the  local  value  at  angle  of  attack  a  .  The  corresponding 
6*  and  6**  can  then  be  evaluated  from  the  definitions 
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The  separation  wedge  based  on  the  equivalent  origin  concept  is  then  utilized  in  carrying 
out  the  separated  flow  analysis. 


f.  Shocks  Imbedded  in  the  Shear  Layer 

In  the  present  model,  a  continuously  spreading  viscous  shear  layer  is  superimposed 
along  the  inviscid  plume  boundary.  At  some  downstream  location,  it  is  quite  possible 
that  the  calculated  internal  boundary  shock  would  be  imbedded  in  the  shear  layer.  This 
occurrence  will  invalidate  the  Chapman-Korst  shear  layer  model  since  normal  pressure 
gradient  and  steep  entropy  rise  are  now  present  within  the  mixing  layer.  A  coupled 
inviscid-viscous  flow  analysis  is  necessary  to  resolve  this  problem.  Complex  mathe¬ 
matical  tasks  are  involved  to  reach  a  solution  and,  to  date,  quantitative  results  have 
not  yet  been  obtained  (Ref.  38). 

5.  PLUME-INDUCED  BOUNDARY  LAYER  SEPARATION 

The  discussion  presented  in  the  foregoing  sections  describes  the  sequence  of  events 
that  take  place  from  the  combustion  chamber  to  the  plume  flow  field  as  well  as  various 
important  jphenomena  pertinent  to  the  LRV  jet  pluming  problem  [flow  regions  (1),  (2), 

(3)  and  (5)  j  .  Because  of  high  jet-to-ambient  pressure  ratios  involved  during  the 
aerocruise  maneuver,  a  relatively  large  equivalent  plume  body  is  frequently  formed 
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which  will,  in  turn,  induce  an  adverse  pressure  gradient  on  the  vehicle  aft-surface. 

If  the  adverse  pressure  gradient  is  sufficiently  strong,  the  boundary  layer  may  not 
have  enough  kinetic  energy  to  overcome  the  retarding  force  and  is  then  forced  to 
separate.  Pressure  elevation  due  to  flow  separation  will  affect  the  aerodynamic  forces 
and  moments  over  the  vehicle  surface.  Serious  stability  problems  may  subsequently 
arise,  which  must  be  adequately  resolved  in  order  to  ensure  satisfactory  flight 
characteristics. 

This  description  shows  that  the  mechanism  of  plume-induced  flow  separation  is  basi¬ 
cally  the  same  as  that  of  separated  flows  upstream  of  a  compression  corner  or  an  aft 
flare,  except  that  mass  exchange  between  the  plume  and  the  air  stream  is  absent  in  the 
solid  problem.  Certain  physical  quantities  between  these  two  types  of  separation  are 
expected  to  be  at  variance.  Semiempirical  formulas  generated  for  the  solid  problem, 
such  as  plateau  pressure  prediction,  may  not  be  applicable  to  the  plume  problem  with¬ 
out  any  modification. 

To  resolve  the  plume-induced  separation  problem,  one  approximate  approach  is  to 
extend  the  Chapman-Korst  shear  layer  model,  originally  developed  for  base  pressure 
predictions  (Refs  26,  27  and  par.  4  of  this  section),  to  the  case  involving  fluid  com¬ 
pression,  mass  diffusion,  and  chemical  reactions.  This  modified  Chapman-Korst 
theory  provides  a  scheme  for  predicting  the  overall  properties  within  the  separated 
flow  region  and  is  merited  with  relative  simplicity. 

For  plume-induced  flow  separation,  the  separation  point  is  also  an  unknown  quantity 
whereas  it  is  usually  prescribed  to  occur  at  the  base  periphery  for  the  base  flow  prob¬ 
lem.  Hence,  an  additional  condition  must  be  Imposed.  In  this  study,  this  condition 
is  provided  by  introducing  a  semiempirical  formula  correlating  the  separation  point 
with  base  pressure,  initial  Mach  number,  and  initial  Reynolds  number.  However, 
development  of  a  versatile  and  satisfactory  formula  for  separation  length  prediction 
still  awaits  the  availability  of  pertinent  data. 

The  modified  Chapman-Korst  model  used  in  LRV  application  is  described  in  the  follow¬ 
ing  section  and  Appendix  I.  Rigorously  speaking,  this  model  is  valid  only  for  two- 
dimensional  or  axisymmetric  flows.  Its  applicability  to  the  three-dimensional  LRV 
plume  problem  in  a  pseudo-two-dimensional  manner  must  be  justified  a  posteriori. 

This  facet  and  the  incorporation  of  the  small  cross -flow  concept  in  the  analysis  are 
discussed  in  par.  6  of  this  section.  The  present  discussion  is  restricted  to  local 
two-dimensional  flow  separation  with  respect  to  a  particular  streamline. 

a.  The  Modified  Chapman-Korst  Model 

In  accordance  with  the  Chapman-Korst  model  (Figs.  1  and  38),  the  plume-induced 
separated  flow  structure  can  be  elucidated  in  a  phenomenological  manner.  When  the 
flow  encounters  a  strong  plume-induced  adverse  pressure,  it  will  leave  the  surface 
and  create  a  dissipative  layer  which,  along  with  the  plume  shear  layer,  will  envelope 
a  region  of  mass  entrainment.  Although  recirculation  occurs  in  this  region,  the 
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Figure  38  Plume-Induced  Separation 
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average  mean  velocity  therein  can  be  regarded  as  very  small  (dead-air).  Hence,  the 
complex  separated  flow  field  is  now  replaced  by  a  comparatively  simple  shear  layer 
envelope.  Implementation  of  mass  and  other  flux  conservations  within  the  shear  layer 
envelope  as  well  as  a  proper  reattachment  criterion  will  then  permit  a  steady-state 
solution  describing  the  average  separated  flow  characteristics. 

The  basic  assumptions  used  in  the  present  analysis  are: 

•  Perfect  gas 

•  Prandtl  and  Schmidt  numbers  are  unity 

•  A  linear  viscosity-temperature  relation 

•  Constant -pres  sure  field 

As  already  discussed  in  the  Introduction  and  par.  4  of  this  section,  as  well  as  in  Ref¬ 
erences  26  and  27,  these  assumptions  are  introduced  for  both  laminar  and  turbulent 
separations.  Review  and  application  of  these  shear  layer  theories  are  also  included 
in  References  39  to  41. 

In  evaluation  of  the  separated  flow  properties,  the  main  flow  quantities  to  be  ascertained 
are  identified  to  be  the  following: 

•  Location  of  separation  point,  x_ 

•  Average  pressure  in  the  plume-induced  separation  region  (base  pressure),  PB 

•  Average  gas  temperature  in  the  separation  region  (base  temperature),  Tfi 

•  Average  species  concentration  in  file  separation  region  (base  species  con¬ 
centration),  Yj,  resulting  from  diffusion  through  the  air  and  the  plume  shear 
layers 

•  The  dividing  (separating)  streamline  location  measured  from  the  shear  layer 
coordinate,  x  (Fig.  38),  defined  j  the  streamline  above  which  the  total 
mass  flux  is  equal  to  the  total  in\  *<cid  mass  flux  immediately  upstream  of 
the  separation  point.  Although  two  dividing  streamlines  are  present,  i.e. , 
one  in  the  stream  shear  layer  and  the  other  in  the  plume  shear  layer,  the 
constant-pressure  condition  permits  a  simple  correlation  between  the  two 
streamlines  (see  Appendix  I). 

•  The  discriminating  (limiting)  streamline  location  measured  from  the  shear 
layer  coordinate,  X  (Fig.  38),  defined  as  the  streamline  which  possesses 
just  sufficient  energy  to  overcome  the  overpressure  in  the  trailing  shock 
system  created  by  the  shear  layer  intersection.  Analogous  to  the  dividing 
streamline  problem,  the  stream  and  plume  discriminating  streamlines  are 
related  in  a  simple  manner. 
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REPRODUCED  FROM  REF.  43 


Figure  39  Reynolds  Number  Effect  on  Plateau  Pressure  Coefficient 


In  view  of  these  six  unknowns,  six  independent  conditions  are  needed  to  obtain  steady- 
state  solutions.  Mathematical  details  on  formulation  of  these  supplementary  conditions 
are  discussed  in  Appendix  I,  and  this  section  places  emphasis  only  on  aspects  that 
attach  particular  physical  significance  to  the  present  LRV  application: 

•  Semiempirical  Relation  Between  X  and  PR.  Inasmuch  as  Pg  represents  an 
average  quantity,  it  can  be  approximated  By  the  average  plateau  pressure  within 
a  separated  flow  region.  A  simple  correlation  formula  for  Pg  can  then  be 
expressed  in  the  widely  used  form  for  laminar  flows  (Ref.  42): 


P_  -  P 
B  o 


2pouo 


(28) 


where  Cpg  is  the  base  pressure  coefficient,  Mq  and  Re0  are  the  Mach  number 
and  Reynolds  number  at  the  separation  point,  and  A  is  an  empirically  deter¬ 
mined  constant.  From  the  various  values  of  A  summarized  in  Figure  39,  a 
constant  value  A  =  1.8  has  been  proven  adequate  for  solid  compression  corner 
flows  and,  from  limited  schliercn  data  (Refs.  2,  6,  and  43),  A  =  G.  5  was 
found  to  be  valid  for  plume-induced  separated  flow.  However,  recent  cone- 
cylinder  pressure  data  with  an  air  plume  at  freestream  Mach  numbers  of  4.5 
and  6.0  (Ref.  25)  have  indicated  that  the  A  value  should  be  much  less  than  6.5. 
Furthermore,  the  NASA-Langley  data  also  show  that  A  is  strongly  dependent 
on  Pex/P0  (jet-to-ambient  pressure  ratio)  but  not  sensitive  to  Mach  number 
for  the  moderate  air  speeds  used  in  the  tests.  Based  on  these  data,  an  empir¬ 
ical  relation  describing  A  is  tentatively  established  and  utilized  in  the  present 
investigation: 


which  is  graphically  represented  in  Figure  40.  The  base  pressure  formula, 
Equation  (28),  essentially  establishes  a  separation  wedge  from  which  impor¬ 
tant  quantities,  such  as  separation  angle,  and  separation  pressure  coefficient, 
can  be  determined  accordingly.  For  turbulent  separation,  a  similar  plateau 
pressure  formula  with  different  values  of  A  and  exponent  for  the  denominators 
can  be  introduced.  However,  no  data  are  available  to  verify  this  conjecture. 

•  Exhaust  Plume  Properties.  For  each  set  of  Xs  and  Pg  values,  the  inviscid 
exhaust  plume  properties  based  on  the  method-of -characteristics  solution  with 
Pg  as  the  local  ambient  pressure  and  under  quiescent  air  conditions  (dead-air) 
can  be  determined  (see  par.  3  of  this  section). 

•  Conservation  Requirements.  Under  steady-state  conditions,  the  total  mass 
Mux ,  th?  mass  flux  for  any  particular  species,  and  the  energy  flux  within  the 
separation  region  must  be  conserved.  In  particular,  formulation  of  the  species 
conservation  requires  that  the  total  species  flux  at  the  separation  point  and 

the  nozzle  exit  equal  the  total  species  flux  at  the  intersection  point  of  the  shear 
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layers  (Fig.  38).  Instead  of  considering  all  the  species  that  exist  in  the 
flow  field,  it  is  more  convenient  to  select  a  particular  species  which  occurs 
in  the  plume  but  not  in  the  air  stream,  such  as  atomic  hydrogen.  By  so  doing, 
the  species  flux  can  be  represented  by  the  exhaust  gas/air  mixture  ratio,  and 
all  other  species  can  be  determined  subsequently  from  chemical  equilibrium 
calculations.  Derivations  of  the  conservation  equations  and  numerical  com¬ 
putation  schemes  are  included  in  Appendix  I. 

•  Thermochemical  Compatibility  Requirement.  The  argument  presented  in 
par.  2,b  of  this  section  has  established  the  validity  of  equilibrium  chemistry 
in  a  dead-air  separation  region.  Thus,  for  a  particular  set  of  base  pressure, 
PB,  and  base  species  concentration,  YB,  a  corresponding  base  temperature, 
TB,  can  be  determined  by  using  the  NAS  A -Lewis  equilibrium  program. 
Furthermore,  the  uniqueness  of  TB  requires  that  the  energy  conservation 
utilizing  the  separated  flow  model  be  compatible  with  that  considered  thermo- 
chemically.  For  H2  -  F2  and  H2  -  O2  systems  with  high-speed  air  flows, 
the  design  charts  presented  in  Figures  2  through  19  can  thus  be  used  as  an 
additional  condition  for  energy  conservation  within  the  separation  region.  In 
application,  TB,  PB,  and  YB  refer  to  the  temperature,  pressure,  and  mixture 
ratio  (CR)  on  these  graphs. 

•  Reattachment  Criterion.  The  reattachment  criterion  requires  that  the  total 

pressure  along  the  discriminating  streamline,  P*  ,  be  equal  to  the  static 

se 

pressure  in  the  trailing  shock  system,  P  ,  so  that  the  flow  can  pass  down¬ 
stream.  The  trailing  shock  system  is  determined  by  the  intersection  charac¬ 
teristics  between  the  separation  wedge  and  the  exhaust  plume  and,  hence,  is 
strongly  dependent  on  the  inviscid  plume  geometry  and  plume  boundary 
properties .  A  simple  two-dimensional  shock  polar  intersection  model  (Ref.  44) 
is  employed  in  this  analysis  to  determine  Pw  as  illustrated  in  Figure  41.  From 
base  pressure,  experiments,  it  has  been  found  that  the  reattachment  criterion  is 
not  always  valid  when  employed  in  conjunction  with  the  Chapman-Korst  model. 
The  difficulty  is  attributable  to  the  over-simplified  model  which  does  not  accu¬ 
rately  describe  the  separated  flow  structure.  To  avert  this  problem,  Nash 
(Ref.  45)  introduced  a  correction  factor  to  validate  the  predicted  values: 


N  =  ^  (30) 

where  P*  is  the  total  pressure  along  the  discriminating  streamline  and  Pw  is 
the  static  pressure  in  the  trailing  shock  system.  The  Nash  factor  is  an  experi¬ 
mentally  determined  parameter;  its  dependence  on  flow  conditions  requires 
further  systematic  study.  Additional  remarks  concerning  this  correction  factor 
are  included  in  Reference  39. 


b.  Computation  Scheme 

With  the  supplementary  conditions  listed  above,  all  the  unknown  parameters  can  be 
determined  from  the  solutions  of  the  conservation  integral  equations.  Simultaneous 


Figure  41  Determination  of  Trailing  Shock  Pressure 
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consideration  of  thermochemistry  and  cross-flow  boundary  layer  would  then  permit 
complete  determination  of  the  plume-induced  flow  field  associated  with  LRV  synergetic 
maneuvers. 

An  iteration  numerical  procedure  has  been  developed  to  ascertain  these  unknown 
parameters ,  tB*  yB\  etc. ,  which  constitute  a  part  of  the  boundary  conditions  for  the 
governing  equations.  The  calculation  procedure  thus  involves  assigning  a  number  of 
values  for  Pb,  Yb,  Tb.  Hse*  etc*  and  evaluating  the  integral  equations  until  all  the 
conservation  and  reattachment  requirements  are  satisfied.  In  performing  the  analysis, 
several  items  are  considered: 

•  To  expedite  the  numerical  iteration,  the  energy  equation  can  be  replaced  by 
the  thermochemical  compatibility  requirement;  then  the  problem  involves  only 
the  relatively  simple  mass  and  species  iterations,  and  Tg  is  correlated  with  Yq 
and  Pb  from  the  thermochemical  design  curves  constructed  independently. 

•  When  the  initial  boundary  layer  thickness  is  finite,  an  apparent  origin  is  defined 
based  on  the  discussion  of  par.  4. 3  of  this  section.  Tlien,  the  corresponding 
separation  wedge  can  be  quite  different  from  that  without  the  initial  disturbance 
(Fig.  37).  TTie  separation  flow  properties  and  the  ensuing  intersection  with 

the  plume  are  determined  in  accord  with  this  apparent  wedge. 

•  Because  of  flat  surfaces  involved  in  LRV  applications  at  least  in  the  sample 
problem  considered  in  this  report,  a  separation  wedge  is  postulated.  It  is 
quite  conceivable,  however,  that  the  separation  region  may  exhibit  conical 
properties,  particularly  for  slender  or  axisymmetric  bodies  (Ref.  43).  Addi¬ 
tional  data  are  needed  to  clarify  this  problem. 

For  mathematical  details  of  the  iteration  procedure,  the  reader  is  again  referred  to 
Appendix  I.  In  this  section,  only  the  computer  flow  diagram  for  the  accelerated  pro¬ 
cedure  (mass  and  species  iterations  together  with  temperature  design  charts)  is  pre¬ 
sented  to  provide  an  overview  of  the  numerical  analysis  (Fig.  42).  Although  a  com¬ 
puter  program  based  on  the  modified  Chapman-Korst  model  has  been  developed,  it  is 
not  an  indispensible  item  since  sufficient  graphical  solutions  have  been  obtained  to  per¬ 
mit  prediction  of  important  parameters,  such  as  separation  distance  and  base  pressure, 
under  a  wide  range  of  flow  conditions.  Application  of  the  pseudo-two-dimensional 
separated  flow  analysis  to  the  LRV  plume  problem  will  be  further  elaborated  in 
Section  m. 

6.  THREE-DIMENSIONAL  BOUNDARY  LAYER  EFFECTS 

The  flow  field  over  an  LRV  surface  is  three-dimensional  in  character  because  of  its 
three-dimensional  configuration  and  possible  pitch,  yaw,  and  roll  motions  during 
maneuvering.  Three -dimensionality  complicates  the  mathematical  analysis  because 
of  the  presence  of  the  cross-flow  component  and  its  derivatives  in  the  governing 
differential  equations.  For  attached  boundary  layer  flows,  the  three-dimensional 
equations  are  amenable  to  similarity  transformation  or  other  analytical  treatments 
under  special  circumstances .  Notable  examples  include  laminar  flow  over  an  infinite 
yawed  cylinder,  the  small  cross  flow  problem,  axisymmetric  spinning-body  flow,  etc. 
However,  their  separated  flow  counterpart  has  not  met  with  comparable  success  be¬ 
cause  of  the  lack  of  understanding  of  the  three-dimensional  separated  flow  mechanism. 
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Figure  42  Computer  Flow  for  Simplified  Scheme 
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Figure  43  Streamline  Coordinate  System  and  Control  Volume 


During  powered  maneuver  of  a  lifting  reentry  vehicle  at  orbital  altitudes,  the  separated 
flow  field  generated  from  the  three-dimensional  attached  flow  by  the  enlarged  exhaust 
plume  would  then  be  three-dimensional  also,  To  predict  the  separated  flow  properties, 
two  questions  naturally  arise: 

•  What  are  the  limitations  of  the  modified  Chapman-Korst  theory  in  application 
to  this  generalized  problem? 

•  Whether  a  simplified  analytical  model  is  feasible,  so  that  the  three-dimensional 
separated  flow  properties  can  be  predicted  approximately. 

In  the  course  of  this  study,  a  cross-flow  boundary  layer  solution  has  been  obtained  which 
not  only  answers  these  two  questions  but  also  allows  rapid  and  realistic  estimates  of 
the  flow  characteristics  associated  with  this  type  of  flow  separation. 

The  basic  concept  underlying  this  cross -flow  theory  is  that  of  perturbation  generation 
for  laminar  boundary  layer  flow  oriented  with  respect  to  the  streamline  coordinate  sys¬ 
tem  (Fig.  43).  Streamline  coordinate  system  means  that  the  x-coordinate  is  coincident 
with  the  external  streamline  so  that  the  cross  flow  at  the  outer  edge  is  zero.  For  a 
large  class  of  boundary  layer  flows,  implementation  of  the  streamline  coordinate  sys¬ 
tem  would  indeed  justify  the  perturbation  approach  for  which  the  cross  flow  can  be  des¬ 
cribed  by  high -order  linearized  solutions.  Inasmuch  as  in  the  current  design  of  lifting 
reentry  vehicles  a  large  portion  of  the  surface  contour  displays  relatively  small  curva¬ 
ture,  e.g.  the  HL-10  surface,  the  small  cross-flow  model  is  expected  to  be  valid. 
Specifically,  the  present  model  evolves  the  following  issues: 

•  It  is  rationalized  that  with  the  streamline  coordinate  system,  the  pseudo-two- 
dimensional  separated  flow  solution  is  valid  along  each  external  streamline. 

•  In  general,  the  external  streamline  is  nit  directly  measurable  but  accessible 
through  surface  oil  flow  observation.  With  the  present  analytic  model,  the 
external  streamlines  can  be  determined  inversely  from  the  small  cross -flow 
solution  together  with  relevant  oil  flow  data.  These  external  streamlines  will 
correspond  to  the  local  streamlines  upstream  of  the  separation  region  dis¬ 
cussed  in  the  preceding  paragraphs,  which  are  needed  for  the  separated  flow 
analysis. 

•  A  frequently  used  technique  to  treat  three-dimensional  boundary  layer  problems 
is  the  strip  theory  for  which  the  vehicle  surface  is  si  bdivided  into  a  number  of 
parallel  strips,  each  exhibiting  a  local  pseudo-two-dimensional  flow  field. 

While  simple  in  principle  and  straight-forward  in  execution,  its  validity  may 

be  limited.  The  small  cross -flow  model  developed  herein  would  permit  a 
critical  evaluation  of  the  strip  theory. 

In  the  following  paragraphs,  therefore,  attention  will  be  directed  to  two  main  topics: 

•  Determination  of  the  external  streamline  from  oil  flow  pattern  and  small 
cross -flow  solution,  including  a  comparison  with  the  strip  theory  results,  and 

•  An  approximate  scheme  to  determine  the  surface  separated  flow  pattern  under 
realistic  three-dimensional  flow  conditions. 
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a.  Small  Cross- Flow  Solution 

On  the  basis  of  the  small  cross -flow  concept  (Refs.  46  to  48),  the  governing  partial 
differential  equations  describing  a  three-dimensional  laminar  boundary  layer  can  be 
reduced  to  a  system  of  ordinary  differential  equations,  the  zero-order  equations 
representing  the  basic  two-dimensional  flow  and  higher -order  equations  describing 
the  cross-flow  effect.  Mathematical  details  of  the  perturbation  development  are  pre¬ 
sented  in  Appendix  II.  It  is  appropriate  to  mention  that  since  all  higher-order  equa¬ 
tions  are  now  linearized,  the  cross-flow  solutions  are  comparatively  easy  to  obtain. 

In  this  analysis,  only  the  first-order  equation  has  been  solved  for  the  cross-flow 
estimate. 

Under  hypersonic  conditions,  the  first-order  equation  can  be  expressed  in  the  following 
form: 

*'0*1"  ^;K-'o2)  ■  °  <3i> 

— * 

where  gj  is  the  cross  flow  velocity  ratio,  p*  is  the  cross-flow  pressure  gradient 
parameter,  (  )'  indicates  differentiation  with  respect  to  the  similarity  y-coordinate , 
and  subscripts  o  and  1  represent  the  zero-order  and  first-order  quantities,  respec¬ 
tively.  Since  the  zero-order  values  are  already  known,  equation  (31)  is,  in  effect, 
a  linear  differential  equation  in  gj.  The  cross-flow  solution  provides  information  for 
flow  distortion  within  a  three-dimensional  boundary  layer.  In  the  present  case,  this 
solution  is  utilized  to  ascertain  the  external  streamlines  with  the  aid  of  appropriate 
oil  flow  data.  The  following  example  demonstrates  this  inverse  cross -flow  application 
and  compares  with  the  predicted  results  using  the  strip  theory. 

Typical  surface  oil  flow  patterns  in  terms  of  angle  of  attack  for  a  blunted  delta  slab 
under  hypersonic  flow  conditions  are  depicted  in  Figure  44  (Ref.  49).  It  is  seen  that 
as  the  angle  of  attack  increases,  the  surface  streamlines  will  bend  outward,  and, 
eventually,  the  geometrical  leading  edge  will  behave  like  a  trailing  edge,  hydrodynam- 
ically.  However,  the  present  concern  is  centered  on  interpretation  of  the  surface 
three-dimensional  flow  behavior  through  the  cross -flow  boundary  layer  theory  rather 
than  evaluation  of  the  angle -of -attack  effect  on  cross-flow  generation  or  the  overall 
aerodynamic  properties.  As  such,  the  zero  angle -of -attack  case  illustrated  in 
Figure  44  can  be  singled  out  for  cross-flow  analysis  without  any  loss  of  generality. 

For  detailed  considerations,  the  zero  angle -of-attack  oil  flow  pattern  at  Mach  6.8  is 
reproduced  in  Figure  45  in  which  three  streamlines  with  azimuthal  angles  of  6,  10, 
and  14  degrees  were  selected  for  investigation.  The  corresponding  P*  values,  where 
(see  Appendix  II) 

f  =  -  (32) 

pl  „  u2  sr  “  2  ej 

C  G 

were  found  to  be  0.02,  0.035,  and  0.05. 
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a  =  20° 


Figure  44 


a  =  25° 


a  =  30° 


a  =  40° 


Surface  Oil  Flow  Pattern  on  the  Compression  Side  of  a 
Blunted  Delta  Slab  at  Mach  6. 8 
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STRIP  THEORY  MODEL 
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Figure  45  Comparison  of  Outer-Edge  Velocities  Based  on  Strip  Theory  and  Small  Cross-Flow  Theory 


The  freestream  flow  was  considered  to  be  of  the  Blasius  type  and  the  cross -flow 
profiles  were  assumed  to  be  approximately  represented  by  those  presented  in  Figure  46. 
Using  the  resultants  of  the  stream  wise  and  lateral  velocity  components  at  n  =  0.  5 
(boundary  layer  thickness  =  9)  to  represent  the  surface  streamlines,  the  inviscid 
streamline  locations  can  be  determined.  The  results  are  depicted  in  Figure  45  and  are 
summarized  below:* 


Streamline 

(I) 

Surface  streamline 
location  (oil 
flow,  deg) 

(ID 

Angle  between  surface 
and  inviscid 
streamline 
(oil  flow,  deg) 

(III) 

Deviation  of  inviscid 
streamline  from 
strip  theory 
result  (deg) 

A 

6 

3.5 

2.5 

B 

10 

6.3 

3.7 

C 

14 

8.7 

5.3 

As  a  clarification  of  the  tabulation  above  and  Figure  45,  prediction  procedure  for 
streamUne  A  is  exemplified  below: 

(1)  Identify  the  surface  oil  flow  angle  for  streamline,  i.e. ,  0'  =  6°  measured 
from  the  center  line,  Column  I  of  the  table. 

(2)  From  the  cross-flow  solution  and  the  oil-flow  streamline,  the  streamwise 
flow  direction  is  found  to  be  0  =  3.5°  measured  from  the  oil-flow  stream¬ 
line  (Column  II).  As  mentioned  above,  n  =  0.5  is  used  in  determining  the 
surface  streamline  components.  The  streamwise  flow  direction  also 
represents  the  local  inviscid  streamline  direction  in  view  of  the  streamline 
coordinate  system  employed  here. 

(3)  Since  the  outer -edge  velocity,  u^,  based  on  the  strip  theory  is  6  degrees 
measured  from  the  center  line,  the  discrepancy  between  the  cross-flow 
method  and  the  strip  theory  is  2. 5  degrees  (Column  HI). 

Limited  results  have  indicated  that  the  discrepancy  resulting  from  the  strip  theory 
model  generally  would  not  induce  any  significant  effect  on  predicted  surface  separated 
flow  patterns  (less  than  10  percent).  This  then  justifies  the  strip  theory  model  for  a 
considerable  latitude  of  separated  flow  applications.  On  the  other  hand,  the  anomalous 
flow  behavior  near  the  swept  leading  edge  or  over  highly  curved  surfaces  could  con¬ 
ceivably  incur  serious  inaccuracy.  However,  an  in-depth  study  of  this  problem  re¬ 
quires  additional  oil  flow  data  which  are  not  presently  available. 


*The  results  are  based  on  a  Blasius  profile  with  zero  heat  transfer.  In  the  aft  portion 
of  the  delta  surface ,  experimental  data  have  shown  that  these  conditions  are  nearly 
attained  (Ref.  49). 
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Figure  46  Typical  Similar  Cross-Flow  Velocity  Profile 
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b.  Surface  Separated  Flow  Pattern 

The  surface  separated  flow  pattern  or  separation  footprint  induced  by  an  LRV  plume 
can  be  determined  by  the  following  simple  scheme: 

•  Divide  the  vehicle  surface  into  a  number  of  parallel  strips  and  consider  each 
strip  individually  (Fig.  47). 

•  Use  the  strip  that  encloses  the  plane  of  symmetry  as  the  reference  strip  and 
perform  plume -induced  separated  flow  predictions  in  accordance  with  the 
theory  discussed  in  par.  5  of  this  section. 

•  Consider  a  particular  strip  away  from  the  plane  of  symmetry  and  define  the 
attached  inviscid  streamline,  Mi,  which  is  parallel  to  the  reference  stream¬ 
line,  M0;  small  cross-flow  correction  may  be  applied  to  obtain  this  compo¬ 
nent.  For  the  simple  configuration  considered  here,  is  equal  to  M0. 

•  Define  the  local  plume  boundary  B.  Because  of  the  dead-air  hypothesis,  a 
uniform  pressure  prevails  within  the  separation  region;  hence,  the  local 
separation  pressure  corresponding  to  boundary  B  is  equal  to  that  for  the  plane 
of  symmetry.  In  other  words,  the  separation  wedges  corresponding  to  dif¬ 
ferent  strips  are  geometrically  similar. 

•  Assume  the  local  trailing  shock  pressure  to  be  equal  to  that  with  respect  to  the 
plane  of  symmetry.  Based  on  this  assumption  and  the  plume  boundary  Mach 
number,  Mj  cos  ^  ,  the  corresponding  separation  wedge  can  be  determined. 
With  this  criterion  and  the  geometrical  similarity  feature,  the  tip  location  and 
the  size  of  the  separation  wedge  can  then  be  predicted. 

•  Repeat  the  process  for  other  strips  and  obtain  the  parabolic  separation  foot¬ 
print  indicated  in  Figure  47.  The  computation  procedure  is  also  illustrated 
in  Figure  48. 
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PLANE  OE  SYMMETRY 


Figure  47  Separation  Wedges  and  Surface  Separation  Contour 


CONSIDER  STRIP  FOR  PLANE 
OF  SYMMETRY 


» 

Figure  48  Computation  Scheme  for  Determination  of  Separation  Footprint 
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#  Section  III 

APPLICATION  TO  A  SIMPLE  CONFIGURATION 

Within  the  framework  of  the  modified  Chapman-Korst  model  development  in  the  present 
study,  the  average  plume-induced  separated  flow  characteristics  during  synergetic 
maneuvering  can  be  predicted.  However,  scarcity  of  relevant  data  has  prevented  full 
endorsement  of  predicted  results  on  a  quantitative  scale.  For  this  reason,  no  effort 
was  made  in  this  phase  of  the  study  to  make  predictions  for  actual  LRV  configurations. 
Instead,  a  simple  lifting  body  configuration  was  thoroughly  investigated  for  the  twofold 
purpose  of  demonstrating  the  analytical  model  and  acquiring  a  qualitative  evaluation  of 
aerodynamic  performance  changes  when  subjected  to  propulsion  effects.  The  results 
will,  of  course,  be  updated  when  additional  data  become  available. 

The  simple  configuration  consists  of  a  delta  planform  with  a  triangular  cross  section, 
as  illustrated  in  Figure  49.  The  delta  flat  plate  is  the  bottom  or  windward  surface 
whereas  the  slope  sides  constitute  the  leeward  surfaces  on  which  the  most  leeside 
ridge  is  at  10  degrees  with  the  delta  surface.  For  the  windward  surface,  all  elemental 
strips  are  assumed  to  be  parallel  to  the  center  line,  as  indicated  in  Figure  45.  For 
the  leeward  surfaces,  it  is  postdated  that  the  elemental  surface  strips  are  parallel  to 
the  most  leeside  ridge  and  are,  therefore,  also  at  10  degrees  with  the  bottom  plate 

•  (Fig.  50).  This  leeside  strip  model  is  consistent  with  that  proposed  in  Reference  50 
and  shown  in  Figure  51.  Although  curved  surface  streamlines  near  the  leading  edge 
are  exhibited  in  Figure  51,  straight  parallel  strips  are  adopted  herein  (Fig.  50)  for 
two  reasons: 

•  The  aim  of  this  study  is  to  assess  the  plume  effect  on  LRV  aerodynamic 
characteristics.  The  lack  of  rigor  for  streamline  determination  near  the 
leading  edge  is  not  expected  to  nullify  the  validity  of  this  assessment,  par¬ 
ticularly  if  the  separation  footprint  is  located  in  the  aft  portion  of  the  surface. 

•  In  predicting  the  average  pressures  over  these  flat  surfaces ,  the  Newtonian 
flow  is  employed,  which  is  compatible  with  the  "straight"  strip  model. 

The  ensuing  paragraphs  encompass  a  detailed  analysis  of  the  plume-induced  separa¬ 
tion  phenomena  associated  with  this  simple  delta-slab  LRV  configuration  as  well  as  a 
critical  evaluation  of  the  present  analytical  model  so  that  a  guideline  for  future  studies 
on  plume-induced  separation  and  related  disciplines  can  be  outlined. 

1 .  PREDICTION  OF  PLUME-INDUCED  FLOW  SEPARATION 

The  freestream  conditions  under  which  the  exhaust  plume  effects  were  studied  are: 

Moq  =  6,  10,  15  and  20 

Altitude  =  140,000,  200,000,  230,000,  250,000,  and  280,000  feet 
«  a=0,  10,  and  20  degrees 
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Figure  49  Plume-Induced  Separation  on  Body  with  Delta  Planform 
and  Triangular  Cross  Section 
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' 


Figure  50  Moment  Arms  from  Leeside  Centroid  to  Point  at  60-Percent  Chord 
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REPRODUCED  FROM  AFFDL-TR-G8-IG1 


SURFACE 

STREAMLINES 


BOUNDARY 

LAYER 


SHOCK  WAVE 


Figure  51  Suggested  Upper  Surface  Flow  Model 
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The  basic  step-by-step  analytical  procedure  to  carry  out  the  separated  flow  prediction 
is  summarized  below: 

•  Define  local  freestream  conditions  for  given  Mach  number,  altitude,  and  angle 
of  attack 

•  Perform  plume-induced  separation  predictions  iteratively,  based  on  methods 
discussed  in  Section  II  and  Appendixes  I  and  n.  The  base  pressure  is  defined 
by  the  plateau  pressure  formula,  Equation  (28).  In  addition  to  the  empirical 
Equation  (29)  to  determine  the  plateau  pressure  coefficient,  the  Nash  factor 
for  the  reattachment  criterion  was  also  determined  empirically.  Data  avail¬ 
able  comprise  a  number  of  schliercn  photographs  by  McGhee  (Ref.  25)  using 
the  same  body  configuration  as  that  selected  in  this  study.  It  was  ascertained 
that  the  Nash  factor  should  be  N  =  0.56  to  obtain  steady-state  separated  flow 
solutions.  This  N  value  was  used  for  a  wide  range  of  Mach  numbers  and  alti¬ 
tudes,  but  because  of  limited  Mach  numbers  (M0#=  4.  5  and  6.0)  used  in  this 
test,  extensive  applicability  for  this  particular  N  value  needs  further  veri¬ 
fication. 

•  Determine  separation  footprints  on  vehicle  surfaces  with  the  aid  of  the  small 
cross  flow  solution,  thus  defining  surface  areas  in  which  pressure  elevation 
occurs 

•  Calculate  pertinent  aerodynamic  parameters  and  estimate  plume  effects  on 
vehicle  performance  and  stability 

a.  Important  Features  of  Problem  Solution 

In  carrying  out  the  analysis,  a  number  of  problems  were  encountered  and  certain  fea¬ 
tures  of  the  problem  solutions  were  exhibited ,  as  follows: 

•  Because  of  the  10-degree  vertex  angle  of  the  body,  the  leeside  surfaces  undergo 
a  10-degree  compression  at  zero  angle  of  attack.  However,  at  10-degree 
angle  of  attack,  the  leeside  surfaces  are  effectively  parallel  to  the  free  stream. 

•  Results  presented  in  Figure  52  show  predicted  dimensionless  separation  length, 
Xs/L,  versus  Mach  number  for  various  Pcx/P0  values  or  local  altitudes.  In 
the  figure,  it  is  clearly  manifested  that  the  separation  region  is  enhanced  by 
an  increasing  Pex/po  but  diminishes  with  increasing  Mach  number.  These 
qualitative  trends  are  consistent  with  the  theories  discussed  in  Section  n. 
Because  of  'local"  altitudes  used,  the  angle-of-attack  effect  is  considered 
implicitly.  The  curves  are  re-plotted  in  Figure  53  in  terms  of  Xs/L  versus 
Pex/P0*  Pertinent  results  are  also  compiled  in  Tables  V,  VI,  and  VII. 

•  It  should  be  mentioned  that  for  Pex/po  =  4030  or  at  280, 000  feet  (Figs.  52 
and  53),  no  steady-state  solution  was  obtained  based  on  N  =  0. 56.  From 
schlieren  observation  (Ref.  25),  however,  a  new  N  value  of  1.5  was  ascer¬ 
tained,  and  the  results  are  depicted  in  Figures  52  and  53.  For  altitudes  above 
280,000  feet  and/or  Mach  numbers  higher  than  20,  the  separation  distance 
was  estimated  by  extrapolation,  as  indicated  in  Figure  53. 

•  For  leeside  problems  involving  high  angles  of  attack,  extrapolation  was  neces¬ 
sary  to  predict  the  plateau  pressure  coefficient  and,  subsequently,  the  base 
pressure.  Furthermore,  a  thick  boundary  layer  frequently  accompanied  the 
leeside  expansion  and  necessitated  additional  extrapolation  in  order  to  reach 

a  solution.  This  then  rendered  the  solution  less  certain  and  open  to  speculation. 
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Table  V  Problem  Solution  Summary,  a  =  0  Degree 
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Subscript  o  refers  to  Local  freestream  conditions  upstream  of  separation  point. 


Table  VI  Problem  Solution  Summary,  a  =  10  Degrees 
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Same  as  a  =  0°  case  if  windward  and  leeward  values  are  interchanged. 
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Note:  20°  compression  on  windward  side/10"  expansion  on  Jeeward  side. 


•  The  predicted  separation  footprints  are  presented  in  Appendix  m.  The  simple 
body  configuration  permits  aerodynamic  predictions  in  a  straightforward 
manner. 


b.  Plume  Effects  on  Aerodynamic  Characteristics 


Without  loss  of  generality,  pertinent  aerodynamic  parameters  are  expressed  in  terms 
of  the  windward  delta  surface  area  and,  for  aerodynamic  moments,  an  additional  refer¬ 
ence  length  of  60  percent  chord  from  the  tip  of  the  delta  surface  (Fig.  50).  However, 
for  the  present  simple  delta-triangle  configuration,  the  force  components  for  the  lee- 
side  or  slope  surfaces  are  dependent  upon  the  geometrical  angles  and  the  pressure 
forces.  In  deriving  a  general  expression  for  the  force  determination,  pertinent  geo¬ 
metrical  parameters,  such  as  illustrated  in  Figure  54,  are  used.  The  main  objective 
here  is  to  ascertain,  N ,  the  unit  vector  normal  to  the  slope  surface  at  any  point  on 
the  surface.  From  Figure  54,  the  side  vectors  can  be  expressed  as 


=  fi  +  ak  =  {  (i  +  tanpk) 
8  =  if  +  hT  -iff*  tanej) 


(33) 

(34) 


ITien,  the  unit  vector  normal  to  the  slope  surface  is  simply  the  cross  product  of  S^and 
S2  as  follows: 

^  S  xS.  ^2  ^ 

N  =  — — —  =  — — —  (T  tanQ  tanP  -  7  tanp  -  k  tane  )  (35) 

ls2xsil  lS2xSll 

where  |S2  xSi|  represents  the  magnitude  of  the  cross  product  and  is  equal  to 

jS^xS^j  =  t2  (tan2e  tan2P  +  tan2p  +  tan20  )1'2  (36) 


If  the  pressure  force  normal  to  the  leeside  surface  is 

Fk  -  apk  \ 


(37) 


where  the  subscript  k  provides  a  general  expression  whether  or  not  jet  pluming  occurs. 
Titus,  normal  pressure  force  with  respect  to  the  windward  reference  surface  = 

-Fk  tanp  (tan2e  tan2p  +  tan2p  +  tan2e  f 1/2  (38) 

which  is  the  j  component,  and  axial  pressure  force  = 

F^  tane  tanP  (tan2e  tan2P  +  tan2P  +  tan2e  )  *^2  (39) 

which  is  the  i  component. 
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Two  limiting  cases  are  also  immediately  obvious: 

Normal  pressure  force  =  !  p  q  I  «  (40) 

("*k  1  r 

In  application  to  the  present  plume  problem,  these  equations  can  be  used  for  both  the 
jet -on  and  jet-off  cases.  From  Figures  49  and  54,  it  is  apparent  that  0  =  p  =  -10°  for 
h  =  a.  The  main  aerodynamic  quantities  are  determined  by  the  following  equations: 

(1)  Jet-off 

Windside  normal  force* 

F  =  (P  -  P„)  A  (41) 

o 

Leeside  normal  force* 


Axial  force 

-1/2 

F^  =  (Pq  -  PK  Ajf  tane  tanp  (tan2e  tan2P  +  tan2P  +  tan2e) 
o  * 

=  0.702  tan  10°  (PQ  -  P^  A ( 

=  0.124(Pq  -  P  Jjf  Af  (43) 

For  completeness,  the  skin  friction  contribution  to  aerodynamic  forces 
should  be  included  in  the  foregoing  equations.  However,  according  to  data 
reported  in  Reference  50  in  which  the  same  body  configuration  was  tested , 
the  skin  friction  effect  on  normal  force  component  was  not  pronounced, 
and  the  frictional  axial  force  was  found  to  be  of  the  same  order  of  magnitude 
as  the  pressure  axial  force.  Consideration  of  the  skin  friction  contribution 
would  not  shed  light  on  understanding  the  plume-induced  separated  flow 

♦The  normal  forces  are  directed  toward  the  respective  surfaces  vectorially. 
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mechanism  but  would  introduce  additional  numerical  integration  and  ambigu¬ 
ity,  the  latter  referring  to  the  problem  of  evaluating  skin  friction  character¬ 
istics  within  a  separation  region  for  which  no  satisfactory  simple  theory  is 
currently  available.  Inasmuch  as  the  present  analysis  is  concerned  with 
determination  of  aerodynamic  force  differentials  and  other  incremental  quan¬ 
tities  arising  from  jet  pluming,  it  is  doubtful  if  omission  of  the  sldn  frictional 
forces  would  affect  the  results  at  all.  Hence,  the  skin  friction  aspect  is 
ignored  in  this  analysis. 

The  corresponding  aerodynamic  coefficients  can  be  expressed  as  follows: 
Normal  force  coefficient 


Axial  force  coefficient 


Pitching  moment  coefficient 


1 


«  C 


X  -  F  X'  +  2F.  Z 
w_ 


(44) 


(46) 


(46) 


where  d  is  the  60  percent  chord  location,  X,  X*  and  Z  are  the  moment  arms 
from  the0 different  centroids  to  the  60  percent  chord  location,  Aw  is  the  bot¬ 
tom  delta  plate  surface  area,  is  the  leeside  slope  surface  area,  and 
q«>  is  the  freestream  dynamic  pressure  parameter  with  respect  to  the  delta 
surface,  i.e. , 


a 


■*( 


2 

P  u 

90  00 


(47) 


According  to  the  symbols  used  herein,  subscript  «  means  the  freestream 
value,  whereas  subscript  o  refers  to  local  freestream  state  upstream  of  the 
separated  flow  region  vrith  shock  pnd  flow  expansion  effects  properly  accounted 
for;  subscripts  w  and  i  are  the  windward  and  leeward  surfaces,  respectively. 
Because  of  symmetry  with  respect  to  the  central  plane,  i.  e. ,  plane  OCB  in 
Figure  54,  all  yaw  and  roll  parameters  are  zero. 

(2)  Jet-on 

By  the  same  token,  the  aerodynamic  quantities  showing  plume-induced  sepa¬ 
ration  effects  can  be  evaluated  by  using  appropriate  values  for  pressure  force 
and  wetted  area. 


4 
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VVindside  normal  force  increment 

F  =  (Pn  -  PJ  \ 
w  '  B  owl. 

Leeside  normal  force  increment 


w 


FN  =  ^PB  "  Po>f  tan2p  +  tan2p  +  tan20) 


-1/2 


=  0.  702  <PB  -  P^  Af 
Axial  force  increment 


t 


9  9  9  2  -l/2 

F  A  =  (P_  -  P  ).  A,  tan  0  tan  p  (tan  0  tan  p  +  tan  p  +  tan  0) 

A  D  O  f  10 

*  “• 124  <PD  -  Po>(  \ 

Incremental  normal  force  coefficient 


iCN  “  ^(Fw  -  Fn) 


Incremental  axial  force  coefficient 


AC  .  =  —  F 


Incremental  pitching  moment  coefficient 


ACM  '  rV  I2  Vf  -  FwXf  +  2FAZf 
80  C 


(45) 


(46) 


(47) 


(48) 


(49) 


(50) 


where  Af  and  Af  are  the  separation  footprint  areas  on  the  windside  and 
leeside  plates,  and  Xf,  and  are  the  moment  arms. 

The  results  are  presented  in  Figures  55  through  70  and  compiled  in  Table  VTII.  Impor¬ 
tant  features  of  these  results  are  as  follows: 

•  The  relatively  small  plumes  for  230,000  and  250,000  feet  altitudes  have 
resulted  in  moderate  changes  in  aerodynamic  properties,  as  expected.  How¬ 
ever,  at  higher  altitudes,  say  280,000  feet,  the  large  plume  would  induce 
more  extensive  flow  separation  and,  therefore,  greater  changes  in  pressure 
coefficients  and  L/D. 
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•  In  general,  exhaust  plume  results  in  a  loss  of  L/D.  The  unexpected  increase 
in  L/D  at  low  angle  of  attack  at  280, 000  feet  shown  in  Figure  70  is  an  excep¬ 
tion,  since  in  the  present  example,  the  separated  flow  region  on  the  windside 
surface  is  much  greater  than  on  the  leeside  surface  at  zero  and  low  angles  of 
attack,  thus  creating  a  strong  overpressure  and  a  higher  L/D.  The  trend  is 
reversed  at  high  angles  of  attack,  as  seen  in  Figure  70.  However,  the  erratic 
behavior  exhibited  for  high-altitude  plume  problems  may,  in  part,  be  attributed 
to  the  extrapolation  procedure  required  in  ascertaining  the  separation  length 
(Fig.  53),  which  could  significantly  affect  the  ensuing  separation  footprint 
predictions. 

•  From  Table  VIII,  it  is  seen  that  all  the  aerodynamic  parameters  are  effective¬ 
ly  independent  of  altitudes  for  the  jet-off  case.  Hence,  the  exhaust  plume 
effects  can  be  demonstrated  by  comparing  Cn«  Ca.  and  Cm  values  for  differ¬ 
ent  jet  pressure  ratios,  where  Pcx/Poo  =  0  refers  to  the  jet-off  case.  These 
convenient  design  curves  arc  illustrated  in  Figures  67,  68  and  69. 

•  The  stability  curves,  Cn  vs  Cm>  illustrated  in  Figures  64,  65  and  66  indicate 
that  plume-induced  flow  separation  affects  the  stability  insignificantly  at 
moderate  altitudes  (230,000  and  250,000  feet).  However,  severe  changes  in 
stability  are  evidenced  a*  high  altitudes  (280,000  feet).  The  results  presented 
in  these  figures  show  the  anticipated  qualitative  trend  and  regions  of  instability 
under  particular  a  and  M  conditions. 

•  Because  of  angle -of-attack  consideration  alone  and  because  of  symmetry  with 
the  vertical  plane  (Fig.  54),  the  yaw  and  roll  effects  do  not  occur  in  this 
sample  problem. 

•  The  design  curves  presented  in  Figures  55  through  69  can  be  extrapolated, 
in  a  limited  fashion,  to  predict  aerodynamic  characteristics  at  higher  .ingle 
of  attack  and/or  higher  Mach  number.  For  example,  at  Mach  number  20, 
250,000-foot  altitude,  and  30-degree  angle  of  attack,  the  plume  effects  arc 
estimated  to  be  as  follows: 

CN  (jet-on)  -  CN  (jet-off)  *  0 

C.  (jet-on)  -  C.  (jet-off)  »  0.002 
A  A 

CM  (jet-on)  -  (jet-off)  *»  -0.01 

Because  of  lack  of  pertinent  data,  extensive  effort  has  not  been  expended  to 
investigate  high  angle  of  attack  and  Mach  number  cases  further.  Also,  the 
validity  of  these  design  curves  for  more  complex  body  shapes  needs  experi¬ 
mental  verification. 
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Figure  61  CN>  C and  CM  vs  a  for  Moo 
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Figure  63  C.T,  C  . ,  and  C„  vs  a  for  M<»  =  15  at  280, 000  Feet 
N  A  M 


Mach  Numbers  at  230, 000  Feet 


Moo  -  6 


pex/p« 


0 

321 

860 

4030 


a 


vs  a  for  Various  Mach  Numbers  and  Jet  Pressure  Ratios 


10 


•9 


I 


Figure  68  vs  a  for  Various  Mach  Numbers  and  Jet  Pressure  Ratios 
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Figure  70  Exhaust  Plume  Effect  on  L/D 
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Table  Vm  Predicted  Aerodynamic  Characteristics  for  Simple  Delta 
Slab  Configuration  With  and  Without  Plume  Effects 
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Remark:  Skin  friction  effects  not  considered 
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2.  CRITICAL  ANALYSIS 

The  main  sources  of  erro^  in  the  present  simplified  analysis  to  predict  plume-induced 
flow  separation  are  itemized  as  follows: 

a.  Plateau  Pressure  Correlation 

Lack  of  relevant  pressure  data  constitutes  the  major  source  of  analytical  uncertainty. 
Deduced  plateau  pressures  from  preliminary  NASA- Langley  cone-cylinder  data  may 
be  quite  erroneous . 

b.  Nash  Factor 

The  Nash  factor  was  introduced  to  reduce  the  discrepancy  between  theoretical  and  ex¬ 
perimental  results.  In  essence,  it  serves  to  rectify  the  over-simplification  of  the 
theorectical  model.  Additional  data,  particularly  at  high  Mach  numbers,  are  needed 
to  establish  whether  this  parameter  is  dependent  upon  the  freestream  flow  conditions. 

c.  Nozzle  Wall  Friction 

Because  of  nozzle  wall  viscous  effects,  the  plume  contour  may  be  much  larger  than  the 
predicted  inviscid  plume;  the  flow  properties  along  the  plume  contour  are  then  modified 
accordingly.  The  plume/separation-wedge  interaction  characteristics  may  be  sub¬ 
stantially  affected  also.  Using  the  case  Moo  =  6,  a  -  0  degree,  and  altitude  =  250, 000 
feet  as  an  example,  it  was  found  that  the  Nash  factor  should  be  0. 3  instead  of  0. 56  If 
the  "l/7th-power"  plume  was  used  (see  Fig.  32).  In  view  of  uncertainty  toward  this 
problem,  a  careful  experimental  program  is  needed  to  reach  a  final  conclusion. 

d.  Chemical  Kinetics  in  Nozzle  and  in  Plume 

The  discrepancy  in  nozzle  performance  between  a  real  nozzle  involving  chemical  kinetics 
and  a  pseudo-ideal  nozzle  with  an  average  constant  Y  of  1. 38  could  induce  considerable 
difference  in  plume  properties,  thus  affecting  the  plume-induced  separated  flow  results. 
By  the  same  token,  the  constant  Y  assumption  for  plume  predictions  beyond  the  nozzle 
exit  could  produce  additional  discrepancy  with  the  real  case.  The  cumulative  effects 
of  the  constant  Y  idealization  may  be  quite  substantial  and,  therefore,  warrant  ample 
attention  in  carrying  out  exhaust  plume  analysis. 

e.  Rarefied  Gas  Effect 

At  high  altitudes  or  under  locally  low-pressure  conditions,  the  rarefied  gas  effects 
should  be  incorporated  in  the  analysis  in  two  aspects:  (a)  the  initial  boundary  layer 
prediction,  and  (b)  hypersonic  viscous-inviscid  interaction  near  the  leading  edge.  Both 
will  influence  the  flow  field  upstream  of  the  separated  flow  region  and  will,  therefore, 
affect  the  predictions  of  the  separation  footprint  and  other  quantities.  However,  at 
very  high  altitudes  where  the  flow  is  of  the  free  molecule  type,  plume-induced  separa¬ 
tion  will  not  likely  be  of  practical  importance  because  of  the  very  small  aerodyanmic 
forces  experienced  by  the  vehicle. 
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f .  Three-Dimensional  Flow  Effect 

Although  it  has  been  argued  in  par.  6  of  this  section  'hat  the  three-dimensional  bound¬ 
ary  layer  effect  is  not  severe  for  the  present  sample  problem,  additional  theoretical/ 
experimental  studies  are  needed  to  reach  an  irrefutable  conclusion.  Aside  from  this, 
two  other  factors  will  affect  the  analytical  prediction: 

(1)  At  angle  of  attack,  complicated  cross-flow  phenomena  may  create  surface 
flow  separation  which  is  not  easily  tractable.  For  purpose  of  demonstration, 
a  surface  separation  curve  for  a  pitched  delta  slab  is  exemplified  in 
Figure  71  (Ref.  37). 

(2)  In  carrying  out  the  present  analytical  procedure,  the  assumption  of  wedge¬ 
like  or  cone-like  separation  will  affect  the  laminar  solution  through  the 
Mangier  and  Dorodnitzyn  transformation  variables  but  will  not  affect  turbu¬ 
lent  flows.  As  seen  in  Appendix  I,  a  deviation  factor  of  Vr~  differentiates 
the  wedge  flow  from  the  conical  flow  where  r  is  a  reference  radial  length. 
For  the  sample  problems  considered  herein,  it  was  found  that  insignificant 
error  would  be  incurred  if  r  is  less  than  3  feet. 

g.  Viscosity  Parameter 

Since  Reynolds  number  is  one  of  the  parameters  used  in  reaching  a  separated  flow 
solution,  an  error  in  dynamic  viscosity,  m,  would  then  lead  to  an  erroneous  solution. 
Again,  from  the  viewpoint  of  numerical  computation,  the  effect  arises  from  the  change 
in  the  transformation  variables  when  different  M  values  are  used.  Similar  to^the  r 
problem  discussed  in  the  preceding  paragraph,  the  square-root  quantity, 
constitutes  the  deviation  factor. 

h.  Shear  Layer  Profiles 

The  idealized  Chapman  and  Korst  profiles  are  subject  to  further  scrutiny,  particularly 
if  the  shear  layer  contains  shocks  and  other  discontinuity  surfaces. 

i .  Equivalent  Origin 

The  hypothesis  that  the  actual  separation  wedge  coincides  with  that  initiated  at  the 
equivalent  origin  must  be  verified  experimentally.  The  equivalent  origin  is  innovated 
to  circumvent  the  mathematical  dilemma  encountered  when  the  initial  boundary  layer 
thickness  is  finite.  The  analytical  inaccuracy  due  to  this  hypothesis  should  be  care¬ 
fully  assessed. 

3.  SIMULATION  PARAMETERS 

Based  on  the  analytical  results  contained  in  this  report  and  the  critical  comments 
summarized  in  the  preceding  section,  a  number  of  correlation  parameters  can  be 
established  to  obtain  similitude  between  actual  LRV  configuration  and  wind  tunnel 
subscale  model.  A  brief  discussion  of  these  parameters  is  presented  below: 
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INCIDENCE  (DEG) 

(•^■Measured  At  Trailing  Edge) 

Figure  71  Position  of  Separation  as  a  Function  of  Incidence 
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a.  Mach  Number 

This  parameter  is  limited  by  the  tunnel  capability,  and  a  complete  range  of  Mach 
number  simulation  is  seldom  possible.  However,  for  a  locally  two-dimensional  plume- 
induced  separation  model  proposed  herein,  the  flat-plate  solution  for  determining 
separation  length  would  be  of  practical  use  for  a  large  number  of  problems.  Under 
thi3  circumstance,  the  Mach  number  simulation  problem  can  be  partially  resolved  by 
angle-of-attack  manipulation  of  the  test  plate.  For  example,  pitching  can  enhance  the 
effective  Mach  number  on  the  leeside  of  the  test  from  which  separation  length  and  other 
parameters  can  be  determined  on  the  basis  of  the  present  theory. 

b.  Reynolds  Number 

Again  tunnel  limitations  would  govern  the  extent  of  Reynolds  number  simulation.  In 
accordance  with  the  present  analysis,  the  Reynolds  number  correlation  can  best  be 
implemented  in  the  following  manner: 

Re  (M2-l)  |a  =  Re  (M2-l)  |s  (51) 

where  Re  refers  to  that  based  on  body  length,  and  subscripts  A  and  S  indicate 
the  actual  and  model  geometry,  respectively.  This  equivalence  is  compatible  with  the 
plateau  pressure  expression  of  equation  (28). 

From  the  discussion  of  the  preceding  section,  viscosity  and  body  length  would  affect 
the  theoretical  prediction  of  separation  length  through  the  Lees-Dorodnitzyn  and 
Mangier  transformation  stated  in  Equation  (1-5).  Hence,  another  Reynolds  number 
parameter  can  be  obtained  by  combining  the  transformation  variables  of  Equation  (1-5) 
as  follows: 


\] l0\“eVe'2' ^ 


s 


(52) 


(53) 


where  Re  is  based  on  body  length  l ,  and  the  longitudinal  distance,  x,  and  the  equivalent 
body  radius,  r,  are  assumed  to  be  linearly  proportional  to  t  Equation  (53)  primarily 
applies  to  the  separation  region  where  an  equivalent  separation  wedge  or  cone  can  be 
defined. 
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c.  Jet-to-Ambient  Pressure  Ratio 

For  fixed  nozzle  geometry  and  V,  the  jet-to-ambient  pressure  ratio,  P  /P  ,  is  the 
main  parameter  that  controls  the  similitude  of  the  plume  boundary  characteristics, 
such  as  the  geometric  contour,  boundary  Mach  number,  etc.  From  the  discussion  and 
non  dimensionalized  plume  contour  curves  presented  in  Section  2,  this  correlation 
parameter  is  clearly  implied. 

d.  Jet  Initial  Turning  Angle 

Again  referring  to  the  discussion  of  Section  2,  the  initial  turning  angle  for  jet  plume 
would  affect  the  initial  plume  contour.  Hence,  a  convenient  correlation  parameter  is 
the  initial  angle 

"l ''iiX <54> 

e.  Jet  Momentum  Coefficient 


In  addition  to  the  two  correlation  parameters  cited  above,  Y  provides  an  additional 
constraint,  and  its  pronounced  effect  on  plume  geometry  has  been  extensively  discussed 
in  Section  2.  As  a  generalization,  a  jet  momentum  coefficient  can  be  defined  to  obtain 
Y  correlation  as  well  as  impingement  force  similitude: 

y  P  A  M2  \  yP  A  M2\ 

ex  ex  ex  ex  l  _  ex  ex  ex  ex|  /ri;v 

- 2  /  “  - - - 5/ 

Yoo  Po  Aq  Moo  /A  P*  A0  Moo/S 

where  Afl  is  a  reference  area.  However,  momentum  ratio  is  not  a  sufficient  constraint 
in  itself  and  the  force  similitude  is  assured  only  if  the  plume  geometric  similitude  is 
obtained. 


f .  The  Kawamura  Lambda  Parameter 


The  Kawamura  lambda  parameter 


y  M  /S 
ex  ex 


(56) 


which  was  derived  in  Reference  51  based  on  the  weak  wave  theory  provides  correlation 
with  wave  reflection  at  the  plume  boundary.  This  similitude  is  primarily  applicable 
to  weak  wave  problems  and  its  validity  for  plumes  involving  large  jet-to-ambient 
pressure  ratio  and  strong  shocks  requires  experimental  verification. 
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g.  Knudson  Number 


As  discussed  in  Par.  e  of  Section  n,  the  present  analytical  prediction  is  significantly 
dependent  on  the  initial  boundary  layer  which  is  affected  by  the  rarefied  gas  phenomena. 
Thus,  a  correlation  based  on  the  Knudson  number  can  be  introduced  in  the  following 
form: 


4  -  £r)s  <57> 

While  the  foregoing  parameters  would  aid  in  establishing  test  simulation,  a  satisfactory 
scaling  law  has  not  been  successful  for  two  phenomena:  chemical  kinetics  in  the  nozzle 
or  plume  flow  field  and  particle  flow  effect  although  the  latter  is  not  relevant  to  the 
candidate  liquid  propellant  systems  discussed  in  this  report. 
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Section  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  modified  Chapman-Korst  theory  developed  during  this  study  provides  a  simplified 
numerical-interpolation  scheme  to  predict  local  plume-induced  separated  flow  prop¬ 
erties  within  a  narrow  region  along  the  chord  of  a  lifting  reentry  vehicle  surface  dur¬ 
ing  aerocruise  maneuver.  By  coupling  with  cross-flow  boundary  layer  analysis,  the 
plume  effects  on  LRV  aerodynamic  characteristics  can  be  evaluated.  Although  sample 
calculations  have  disclosed  results  consistent  with  anticipated  trends  and  in  agree¬ 
ment  with  limited  data,  semiempiricism  inherent  in  the  analysis  has  precluded  a  wide 
applicability  of  this  method  due  to  the  lack  of  sufficient  data.  Hence,  in  any  future 
study  the  experimental  aspect  must  receive  the  immediate  attention  after  which  further 
theoretical  improvement  can  be  instituted. 

On  the  basis  of  the  discussions  concerning  theory  development  and  a  critical  assess¬ 
ment  of  the  analytical  model  in  the  preceding  sections,  it  is  recommended  that  the 
following  tasks  constitute  the  bulwark  of  any  future  study  on  plume-induced  flow 
separation: 

a.  Experimental  Phase 

(1)  Plateau  Pressure.  Perform  systematic  surface  pressure  measurements  for 
a  wide  variety  of  freestream  flow  conditions  so  as  to  permit  development  of 
an  adequate  plateau  pressure  (or  average  separation  pressure)  formula. 

(2)  Visual  Observation.  Employ  the  oil  flow  method  to  determine  surface  separa 
tion  pattern;  use  schlieren  and/or  shadowgraph  techniques  to  provide  further 
observation  of  flow  separation. 

(3)  Temperature  Measurement.  Provide  information  for  evaluation  of  temper- 
ature  and  surface  heat  transfer  characteristics  within  the  separation  region. 

(4)  Exhaust  Plume  Properties.  Compare  predicted  and  observed  exhaust  plume 
configurations  in  order  to  evaluate  the  validity  of  the  simple  constant  y  model 
and  to  gain  confidence  in  plume  and  nozzle  flow,  predictions  involving  flow 
separation  upstream  of  base  region.  Different  propellant  systems  (thus  vary¬ 
ing  the  combustion  temperature,  the  specific-heat  ratio,  etc.)  should  be  used 
to  gain  better  insight. 

(5)  Configuration  Effect.  Use  different  configurations,  such  as  those  with  fins 
and  flaps,  to  accrue  knowledge  on  three-dimensional  flow  phenomena  and  to 
permit  generalization  of  the  theory. 

(6)  Simulation  Parameters.  To  aid  in  the  experimental  study,  a  set  of  simula¬ 
tion  parameters  can  be  defined  to  carry  out  a  meaningful  wind-tunnel  test 
program.  These  parameters  are  discussed  in  the  preceding  section. 
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b.  Theoretical  Phase 


(1)  Nash  Factor.  Evaluate  this  reattachment  correction  factor  in  terms  of  free- 
stream  flow  variables  so  that  further  mathematical  refinement  of  the  sepa¬ 
rated  flow  analysis  need  not  be  pursued.  Introduction  of  the  empirically  deter¬ 
mined  Nash  factor  permits  approximate  description  of  the  plume-induced 
separation  phenomena  by  means  of  the  pseudo-two-dimensional  model.  For 
example,  Appendix  I  shows  that  the  hypothesis  of  Nash  factor  of  0.56  is 
necessary  to  solve  the  delta  slab  problem.  The  non-unity  Nash  factor  repre¬ 
sents  an  overall  correction  of  cross-flow  mass  diffusion  and  other  complex 
flow  phenomena  that  characterize  the  delta  slab  flow  separation. 

(2)  Rarefied  Gas  and  Hypersonic  Flow  Effects.  Perform  a  systematic  analysis 
on  the  phenomena  of  boundary  layer  growth ,  its  interaction  with  the  inviscki 
stream,  and  the  ensuing  effects  on  plume-induced  flow  separation  under 
hypersonic  and/or  rarefied  gas  conditions. 

(3)  Vortex  Generation.  Determine  the  leading-edge  vortex  generation  effect  on 
plume-induced  flow  separation.  The  LRV  vortex  problem  is  three-dimensional 
in  character,  since  it  is  strongly  dependent  on  the  angle  of  attack  (also  yaw 

or  roll  motion)  and  the  leading  edge  geometry.  Attention  should  be  directed 
to  engineering  prediction,  such  as  estimation  of  separation  wedge  angle  and 
size  upon  encounter  of  vortex  sheets  in  the  flow  field. 
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Appendix  I 

PLUME-INDUCED  BOUNDARY  LAYER  SEPARATION* 

By  M.  C.  Fong 

The  basic  assumptions  for  the  modified  Chapman-Korst  model  to  predict  plume-induced 
separated  flow  characteristics  are  recapitulated  as  follows: 

•  The  inviscid  flow  properties  upstream  of  the  separation  point  are  known  a  priori. 

•  The  flow  is  two-dimensional  or  axisymmetric.  On  a  lifting  reentry  vehicle  sur¬ 
face,  the  flow  is  considered  to  be  locally  two-dimensional  with  respect  to  a 
particular  streamline  (pseudo-tWo-dimensional). 

•  The  flow  separation  is  either  pure  laminar  (Chapman)  or  fully  developed  turbu¬ 
lent  (Korst). 

•  The  average  flow  velocity  within  the  separation  region  is  very  low  (dead  air) 
and  the  average  pressure  therein  can  be  represented  by  the  plateau  pressure. 

•  The  Prandtl  number  and  the  Schmidt  number  are  unity,  thereby  rendering  the 
generalized  Crocco  relation  valid. 

•  In  the  dead-air  region,  an  average  temperature  and  a  fuel-to-air  ratio  can  be 
defined. 


♦The  basic  analytical  model  was  developed  during  a  Navy-sponsored  study  under  Con¬ 
tract  N0003066C0186  but  modified  and  improved  during  the  present  study  as  well  as 
during  preparation  of  the  technical  proposal  for  this  study  contract.  This  appendix 
represents  an  improved  version  of  Reference  41, 
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As  described  in  considerable  detail  in  the  text,  the  unknown  parameters  associated 
with  plume-induced  flow  separation  are: 

•  Separation  location,  Xg 

•  Base  pressure,  Pg 

•  Base  temperature,  Tg 

•  Base  species  concentration,  Yg 

•  Dividing  streamline  location,  17. 

Je 

•  Discriminating  streamline  location,  17 

s 

e 

lrom  which  six  conditions  are  required  to  ascertain  these  unknown  functions.  As  pro¬ 
posed  herein,  the  required  conditions  are  provided  by  a  set  of  integral  equations  rep¬ 
resenting  conservation  of  fluxes  within  the  separation  region,  supplemented  by  an 
appropriate  plateau  pressure  correlation  formula  and  a  reattachment  criterion,  viz. , 

•  Semiempirical  formula  correlating  P_  with  X 

D  S 

•  Dividing  streamline  determination 

•  Conservation  of  total  mass  flux 

•  Conservation  of  energy  flux 

•  Conservation  of  mass  flux  for  individual  species 

•  Reattachment  criterion  -  Comparison  of  the  total  pressure  along  the  discrim¬ 
inating  streamline  with  the  static  pressure  in  the  trailing  shock  system. 

This  appendix  discusses  these  integral  equations  and  supplementary  conditions  with 
emphasis  on  detailed  mathematical  formulation  and  calculated  procedure. 


1.  ANALYTICAL  PROCEDURE. 


Before  step-by-step  calculation  procedure  is  discussed,  several  useful  concepts  and 
definitions  are  summarized  here  for  clarity  and  convenience  of  discussion: 

•  A  dividing  (separating)  streamline  is  the  streamline  which  separates  the  inte¬ 
grated  mass  flux  upstream  of  the  separation  point  from  the  mass  entrained  in 
the  separation  region  (Fig.  1-1). 


•  A  discriminating  (limiting)  streamline  is  the  streamline  which  contains  just 
sufficient  kinetic  energy  to  overcome  the  trailing  shock  pressure  and  enters 
the  compression  zone. 


•  If  the  initial  boundary  layer  thickness  is  finite,  an  equivalent  origin  based  on 
the  "origin  shift"  model  for  re-establishing  self-similar  shear  layer  is  used. 
At  the  equivalent  origin,  the  mass  and  momentum  fluxes  are  the  same  as  that 
at  the  actual  separation  point  (Refs.  35  and  36).  The  origin  shift  coordinates 
are  found  to  be 


xo 


A**  P  U 
<5 _ ^e _ e 

26  Me 


(laminar) 


yo  =  6*  +  6 ** 


(I-D 


X  *  O  6 

o 


77  a 


<p(  1  -9)  dr? 


(turbulent) 


=  6*  +  6** 


d-2) 


where  6  *  is  the  displacement  thickness,  6**  is  the  momentum  thickness,  <p  is  the 
turbulent  velocity  profile  and  v  =  cry/x  is  a  transformed  y  coordinate. 
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lh 


Figure  1-1  Plume -Induced  Separation 


« 


With  the  assumptions  listed  above,  the  complex  flow  field  can  be  considerably  simplified. 
A  further  simplification  in  this  analysis  requires  that  the  separated  flow  field  be  either 
pure  laminar  or  fully  developed  turbulent.  Specifically,  the  shear  layer  on  the  plume 
side  is  assumed  to  be  always  fully  developed  turbulent;  on  the  other  hand,  because  of 
high-altitude  cases  under  consideration,  the  streamside  shear  layer  is  regarded  as 
invariably  laminar. 

For  laminar  separation,  the  governing  momentum  differential  equation  for  the  shear 
layer  can  be  represented  by  a  Blasius-type  equation: 


where  x  and  y  are  oriented  along  the  separation  wedge ,  and  the  boundary  conditions 
are 


u  =  u  for  y  =  « 
e 

and  u  =  v  =  0  for  y  =  •« 


d-4) 


The  Lees-Dorodnitzyn  and  Mangier  transformations  are  employed  to  modify  the  govern¬ 
ing  differential  equations,  i.e. ,  * 


♦The  symbol  £  is  introduced  here  to  avoid  confusion  with  the  turbulent  transformation 
variable,  17  .  In  the  text,  however,  17  is  used  for  both  laminar  and  turbulent  cases. 
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Also,  a  similarity  function  is  defined  as  follows: 


V  (t>  =  <P  -  -jj-  (1-6) 

e 

where  c  =  0  for  two  dimensional  flow  and  c  =  1  for  axisymmetric  flow.  The  re¬ 
sulting  equation  is  of  the  form 


f"  +  ffti  =  o  (1-7) 

for  which  the  solution  is  presented  in  Figure  34  and  expressed  in  numerical  form  in 
Equation  (15)  in  Section  n. 

For  turbulent  separation,  the  momentum  equation  is  linearized  in  accordance  with  the 
Korst  model  (Ref.  27 )  to  become 


u 


8u 

5* 


d-8) 


where  E  is  an  eddy  viscosity  parameter  and  boundary  conditions  similar  to  that  of 
Equation  (1-4)  can  be  established. 


For  zero  initial  boundary  layer  thickness,  an  error  function  solution  will  result 
(Equation  16),  i.e. , 


'(?)!  «-> 

Following  the  calculation  procedures  discussed  in  References  40  and  41 ,  the  plume- 
induced  separated  flow  characteristics  are  determined  by  the  following  steps: 

I 

1.  A  separation  point  is  assumed.  From  given  initial  flow  conditions  and  by 
invoking  the  compression  surface  analogy,  the  plateau  pressure  within  the 


127 


plume -induced  separation  region  can  be  obtained  from  Equation  28  in  the  main 
text  (laminar  separation  only),  i.e. , 


where  the  appropriate  "A"  value,  though  not  yet  firmly  established,  is  pre¬ 
sented  in  Equation  (21).  Equation  (I— 10)  correlates  Pp  with  the  Mach  number 
and  Reynolds  number  at  the  separation  point.  Because  of  the  l/4th  power  for 
the  exponent,  is  generally  not  sensitive  to  the  different  assumed  values 
of  Xs. 

2.  Based  on  the  plateau  (base)  pressure,  a  separation  wedge  emanating  from  the 
equivalent  origin  can  be  defined  from  which  the  local  outer-edge  flow  properties 
are  determined  (Fig.  1-1).  For  example,  the  wedge  angle,  the  separation 
shock,  etc.  can  be  evaluated.  Tentatively,  a  wedge  model  is  used;  modifica¬ 
tions  will  be  made,  should  future  data  disclose  contrary  evidence. 

3.  The  exhaust  plume  contour  and  pertinent  flow  properties  are  determined  based 
on  the  plateau  pressure  (effective  ambient  pressure)  and  the  MOC  solution. 

The  dead-air  hypothesis  permits  the  plume  definition  under  quiescent  flow 
conditions. 

4.  From  the  interaction  between  the  separation  wedge  and  the  exhaust  plume 
contour,  a  trailing  shock  system  is  established.  The  trailing  shock  charac¬ 
teristics  can  be  ascertained  by  means  of  the  usual  compressible  stream  inter¬ 
action  considerations.  Although  viscous  shear  layers  are  assumed  to  develop 
along  both  the  separation  wedge  and  the  plume  contour,  the  intersection  phe¬ 
nomenon  refers  to  the  inviscid  stream  interaction  only. 

I 

5.  Several  values  of  the  average  base  temperature,  Tp,  and  overall  fuel-air 
mixture  ratios  Yn,  in  the  separation  region  are  next  assumed. 
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6.  In  accordance  with  the  generalized  Crocco  relation  for  Pr  =  1  and  Sc  =  1, 
the  total  enthalpy  and  species  concentrations  can  be  correlated  with  the  velocity 
profile  as  follows: 


Stream  Side  (0  <  Y  <  Yg) 


Plume  Side  (Y„<  Y<  1) 

-  D 

H  ■  (H)  -  hb)  ij  +  hb 


(I’ll) 


d-12) 


(1-13) 


(1-14) 


where  the  fuel-air  mixture  ratio,  Y,  has  the  same  value  as  the  hydrogen  ele¬ 
mental  mass  fraction  for  the  present  problem,  i.e. ,  Yj  =  1  and  Yfi  =  0, 


♦ 
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7.  A  number  of  gas  properties  needed  for  the  iteration  process  are  defined  as 
follows: 


c 

P 


(i  -  Y)  c 

*e 


(1-15) 


or 


CD  =  CD  +  (c  "  cn  )  (Stream  Side) 
F  hB  *e  e 


(I-  15a) 


c  =  c  +  (c  -c  )—  (Plume  Side)  (I-  15b) 

P  PB  Pj  PB  Uj 


i^I) 

m  ’ 


(1-16) 


where  equivalence  of  Equations  (1-15),  (I-15a),  and  (I-  15b)  can  be  easily  verified 
with  the  aid  of  Equations  (1-12)  and  (1-14). 
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(1-17) 


♦ 


(Stream  side) 


(Plume  side) 


(1-18) 


where  cpe/cp  and  cpj/cp  are  obtained  from  Equations  (I-15a)  and  (I-15b),  and 
Cr  is  the  Crocco  number  defined  as 


(1-19) 


8.  Several  values  of  discriminating  streamlines  are  then  assumed.  The  Mach 
number  along  each  external  discriminating  streamline  is  computed  from 


131 


Similarly,  on  the  plume  side,  we  have 


M 

J 


and 


(1-22) 


(1-23) 


Thus,  for  each  value  of  assumed  external  discriminating  streamline,  there  is  a 
corresponding  plume  discriminating  streamline  where 


(1-24) 


Because  of  Equation  (1-24),  only  one  discriminating  streamline,  say  the  stream- 
side  one,  needs  to  be  Introduced  explicitly. 


9.  To  express  the  definition  of  dividing  streamline  mathematically,  the  total  mass 
and  momentum  balances  between  the  initial  point  and  any  downstream  location 
(for  zero  initial  boundary  layer  thickness)  are: 


p  u  y 
'e  e  Jo 


d-25) 
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(1-26) 


y 


2 

p  u  y 
re  e  Jo 


dy 


from  which  the  following  equation  results 


(1-27) 


Now  consider  an  arbitrary  streamline  K,  it  is  obtained  that 


yk 


•/ 

y\ 


pudy  =  pe  ue 


yj 


dy 


(1-28) 


Hence,  along  the  dividing  streamline,  y^  ~  y^  and  =  0,  Equation  (1-28) 
becomes 


4 


r 
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Pu 

u 

ne  e 


dy 


(1-29) 
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Equation  (1-29)  is  applicable  to  both  laminar  and  turbulent  flows.  For  a  given 
problem  where  the  velocity  profile  is  known  and  proper  coordinate  transforma¬ 
tion  is  implemented,  the  dividing  streamline  location  jr  can  be  ascertained 
uniquely: 


Laminar 


F  U)  dt 


oo 

J  f  U)  [1  -  V  (?)]  df 

-OO 


Turbulent 


(1-30) 


in  which  the  denominator  contains  the  transformed  density  ratio  as  follows: 


(1-32) 


In  deriving  Equation  (1-32),  the  definition  of  Cr  (Eq.  (1-19))  and  the  following 
supplementary  expressions  are  used: 


V-  1 
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1-C' 


1 


1-C 
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r 


and 
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+ 
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now 


from  which  Ta/T  can  be  expressed  in  terms  of  T^/Tj  as  follows: 


Ta 

T 


RaTta 

RT. 


-1 


hence 


RT 
a  a 

RT 


RaTta 

RT. 


(1-Cr  ) 
a 

RaTta  (y-l\ 
RTt  [y^ 


which  leads  to  Equation  (1-32)  by  substituting 


RT. 


A  = 


RaTta 


10.  For  conservation  of  mass  flux  per  unit  width,  it  is  asserted  that  the  mass  flux 
above  the  discriminating  streamline  can  enter  the  trailing  shock  system.  For 
example,  on  the  stream  side,  we  obtain 


I 


oc 


Pu  dy  =  re 


p  u  dy 


+ 


Pu  dy 


(1-33) 
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Similarly,  on  the  plume  side 


(1-34) 


For  a  steady-state  solution,  mass  flux  conservation  requires  that  the  sum  of 

m  and  m.  be  equal  to  the  total  flux  entering  the  separation  region  plus  mass 
®  J 

bleed  if  it  exists.  In  other  words, 


P 

+  mj  ■  mje  +  mjj  +  mB  *  D  /  Pu  dy  +  mje  +  m 


JJ 


which  results  in 


(  f 

I  Pudy+/  Pu 

ys>  yBl 


dy  =  m 


B 


(1-35) 


This  equation  exhibits  a  mass  corridor  as  if  the  flow  were  entering  in  and  exit¬ 
ing  from  a  corridor  bounded  by  the  dividing  streamline  and  the  discriminating 
streamLie  (Fig.  1-1).  As  a  corollary,  it  is  seen  that  for  the  ordinary  base 
flow  problem  (no  jet  plume  and  no  base  bleed) , 


(1-36) 


indicating  that  the  dividing  streamline  and  the  discriminating  streamline  are 
coincident. 
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In  terms  of  transformed  coordinates,  the  mass  flux  integrals  can  be 
expressed  as: 

•  For  the  laminar  stream  shear  layer 


(1-29) 


d-37) 


(1-38) 


(1-39) 
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where  A  =  (R  T./R  T»  ),  <s>  =  (u/u  ),  subscript  a  indicates  some  references 

l  &  g  6 

state*  and 


‘i  <VV  = 


rj.  .  . 

T  -  c0 

2/y  -  1  \ 

•V  (H 


(1-40) 


The  total  mass  conservation  then  reads 


£  *  -  f  Fl<V{Se>  V  -  *B  ■  0 


(1-41) 


11.  The  total  energy  flux  of  the  separation  region  consists  of  convention  and 
conduction  contributions:  In  the  laminar  case. 


(1-42) 


*The  same  Korst  velocity  profile  and  spreading  parameter  are  used  for  both  two- 
dimensional  and  axisymmetric  cases. 
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Or 


Je  =  & 

rE  r< 


_<He  -  HB»  r2  (v  «..)  +  »B  F!  (£V  £se) 

+  <He-HB»  F1  (£V-)  -  <He  *  V  F2(£l0"") 


(1-43) 


where  the  Crocco  relation  is  used  in  the  equations,  and 


F2«.-*b> 
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Its  turbulent  counterpart  is  found  to  be 


•  • 


E.  =  E  +  E  . 
t^  conv  cond 
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(1-44) 


(1-45) 


(1-46) 
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and 


*3  (V  nb* 


(■  - O 


dr? 


The  energy  flux  conservation  thus,  is  expressible  as  follows: 


fiHBFl(tJe't0  +  <He-HB,|F2(V£so) 


+  F1<V">  •  F2  V"> 


+  J<rj  1  jHB  [l(v  "sj  +  (H)  "  »B>  %) 


(1-47) 


+  Ii(Tjjj*#o) "  w*  lr  qb  =  0 


(1-48) 


where  Qg  is  a  base  heat  transfer  parameter. 


12.  Similarly,  the  species  conservation  for  laminar  and  turbulent  shear  layers 
are  determined  as  follows: 


;?-/  puy*-5f  fYBri(-( 

ysj  L  ' 
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YB>  F2 
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(1-49) 


140 


and 


Y 

J 


Y  dy 


“M-v) 


<j  -  v  y-v 


) 


(1-50) 


Inasmuch  as  the  propellant  gas/air  ratio  is  used  to  represent  Y,  we  have 

Y„  =  0  and  Y.  =  1.  It  follows  that  since  the  species  flux  at  the  initial  sta- 
c  J 

tions  must  be  equal  to  that  flowing  downstream ,  the  following  relation  must 
be  true 


*• 

Y 

e 


In  terms  of  transformed  coordinates ,  we  have 


(1-51) 


£  Y=  —  Y 
rf  B 

Fi('V”) 

-r*M 

p.  U.  X.  . 

a  j  ‘‘  (V  "si) 

*  j 

/  \ 

1  +  «*  -  YB>  ‘2  (-•  %J 

0  (1-52) 


In  summary,  equations  (1-41),  (1-48)  and  (1-52)  are  the  conservation  equations 
employed  in  this  analysis  for  numerical  computations.  As  stated  earlier, 
these  equations  require  laminar  stream  shear  layer  and  turbulent  plume  shear 
layer.  If  both  the  stream  side  and  the  plume  side  are  fully  developed  turbulent, 
the  following  system  of  integral  equations  are  obtained: 
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Mass  Flux 


13.  For  a  given  Pg  and  several  assumed  values  of  Tg  (or  Hg),  Yg,  and  rjg  , 
the  ci  culated  results  of  mass  flow,  energy  flux,  and  species  concentration  can 
be  presented  graphically,  as  depicted  in  Figure  1-2.  Cross-plotting  the  inter¬ 
section  points  will  permit  interpolation  of  these  results  to  ascertain  the  desired 
values  of  Tn ,  Y_ ,  and  n  for  which  the  conservation  requirements , 

yiifa  =0,  5^  E  =  0,  and^Y  =  0  are  fulfilled;  that  is.  Equations  (1-41),  (1-48), 
and  (1-52)  or  Equations  (1-53)  through  (1-55)  are  satisfied. 
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B 
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Figure  1-2  Interpolation  Scheme  for  Determination  of  TB  and  Yg 
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Since  the  gas  mixture  in  the  separation  region  is  chemically  reacting,  the 
separated  flow  parameters,  Tg,  Yg,  etc. ,  predicted  by  the  foregoing  pro¬ 
cedure,  must  be  compatible  with  the  thermochemical  results.  For  example, 
if  chemical  equilibrium  prevails,  then  the  predicted  Tg  and  Yg  should, 
within  the  accuracy  of  the  numerical  method,  agree  with  that  by  means  of  the 
NASA- Lewis  equilibrium  program.  Under  certain  conditions,  such  as  pro¬ 
nounced  chemical  kinetics  effects,  additional  iteration  parameters,  e.g. , 

Cp  or  y ,  may  have  to  be  introduced  in  order  to  determine  all  the  pertinent 
flow  quantities.  However,  this  facet  has  not  been  sufficiently  explored  to 
permit  any  irrevocable  conclusions. 

14.  Following  the  determination  of  Tg,  Yg,  and  Vs  ,  the  stagnation  pressure 

along  the  discriminating  streamline,  Pt  ,  is  compared  with  the  wake  pressure 

s 

in  the  trailing  shock  system,  Pw.  The  method  for  determining  P^  is  discussed 
in  Section  II  and  an  interpolation  scheme  is  graphically  represented  in  Figure  41. 
From  recent  studies  (Refs.  39  and  45),  it  has  been  concluded  that  a  semiempirically 
obtained  Nash  factor  is  generally  necessary  so  as  to  establish  a  good  pressure 
correlation.  The  Nash  factor  is  defined  as 


N 


(1-56) 


which  ranges  from  0.3  to  1.0  for  most  practical  applications.  For  a  delta 
plate  with  triangular  cross  section  discussed  in  the  text,  N  =  0.56  was 
found  to  be  valid,  at  least  for  moderate  Mach  numbers  and  altitudes.  Essen¬ 
tially,  the  Nash  factor  is  devised  to  improve  the  oversimplified  model  in  an 
attempt  to  describe  the  very  complex  separated  flow  mechanism,  and  it  pro¬ 
vides  a  supplemental  condition  for  an  engineering  estimate  of  plume-induced 
separated  flow  properties.  Additional  experimental  information  is  needed  to 
define  the  Nash  factor  adequately  and  to  establish  the  influence  of  flow  condi¬ 
tions  on  this  parameter. 
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If  Ptg  does  not  correlate  with  Pw,  then  a  new  base  pressure,  Pg,  will 
have  to  be  assigned.  The  iteration  process  from  Step  (1)  to  Step  (14)  will  then 
be  repeated.  Pressure  compatibility  thus  provides  the  final  step  for  predic¬ 
tion  of  the  flow  characteristics  associated  with  plume-induced  separation. 

2.  DISCUSSION. 

A  computer  program  for  UNIVAC  1108  has  been  developed  to  carry  out  the  iteration- 
interpolation  process  for  the  plume-induced  separated  flow  solution  discussed  above. 

A  flow  diagram  for  the  computer  program  is  shown  in  Figure  1-3,  and  a  more  expedi¬ 
ent  procedure  using  the  thermochemical  compatibility  requirement  in  place  of  the 
energy  iteration  is  presented  in  Figure  42.  It  should  be  emphasized  that  this  com¬ 
puter  program  was  used  in  generating  proper  design  curves  so  that  the  reader  can 
estimate  the  plume-induced  separation  effects  on  LRV  aerodynamic  performance  and 
stability  semi -graphically.  These  curves  refer  to  the  separation-distance  curves  de¬ 
picted  in  Figures  52  and  53. 

Nevertheless,  in  clarifying  the  analytical  scheme,  two  sample  problems  are  pre¬ 
sented  in  this  appendix.  The  first  example  is  that  of  the  solution  for  a  cone -cylinder - 
flare  problem  is  presented  here  to  demonstrate  the  present  analytical  model.  The 
problem  is  in  accord  with  a  particular  NASA- Langley  air  plume  test  with  =  4.65, 

Re  =  346500,  and  zero  angle  of  attack  (Ref.  2).  The  plume  boundary  was  traced  from 
the  schlieren  photograph;  base  heating  and  diffusion  effects  were  neglected.  On  the 
basis  of  unit  Nash  factor,  the  predicted  results  in  terms  of  separation  point  locations 
arc  depicted  in  Figure  1-4.  Good  agreement  between  predicted  and  measured  values 
is  clearly  indicated. 

In  application  to  problems  involving  thermochemical  considerations,  the  alternate 
scheme  based  on  thermochemicai.  compatibility  (Fig.  42)  is  employed  to  avoid  tedious 
cross -plotting.  Here,  appropriate  equilibrium  thermochemicai  design  curves  are 
constructed  independent  of  the  separated  flow  analysis  but  utilized  in  conjunction  with 
the  iteration  scheme  through  the  uniqueness  of  Tg  ana  Yg  under  equilibrium  conditions. 
For  further  clarification  another  sample  problem  is  included  in  this  appendix  to 
exemplify  the  alternate  computation  procedure.  Selected  for  demonstration  is  the 
simple  configuration  discussed  in  the  main  text  (delta  planform  with  triangular  cross 
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Figure  1-3  Computer  Flow 
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section  with  a  length  of  45  feet  for  the  bottom  delta  plate),  and  the  solution  is  only  for 
that  of  a  thin  strip  containing  the  plane  of  symmetry  wherein  the  flow  is  assumed  to  be 
locally  two-dimensional.  Other  requirements  are: 

•  Windward  Surface  (bottom  plate) 

•  Altitude  =  250,000  feet,  =  0.0486  psf 

•  Mach  Number  =  10 

•  Zero  angle  of  attack;  P  =  P 

O  w 

•  Hydrogen-fluorine  jet  plume  with  P  =  500  psi,  all  other  properties  likewise 
being  consistent  with  that  discussetf  in  the  Section  2. 3. 1. 

The  step-by-step  procedure  with  the  aid  of  proper  illustrations  is  as  follows: 

(a)  Select  three  separation  points;  X  =  20,  30,  and  40  ft  from  the  tip 

O 

(b)  The  resulting  PB  based  on  Equations  (28)  and  (29)  in  par.  5a  of  Section  II 
was  found  to  be  0. 185  psf  for  xg  =  30  ft.  For  simplicity,  this  Pg  value 
was  used  for  xg  =  20  and  40  ft  also.  Because  of  the  l/4th-power  for  the 
denominator,  no  significant  error  would  be  incurred  by  this  simplification. 

(c)  Determine  exhaust  plume  contour  and  other  boundary  properties  using  this 
Pg  as  the  effective  ambient  pressure,  and  zero  freestream  velocity  compati¬ 
ble  with  the  dead-air  hypothesis. 

(d)  Assign  three  base  temperatures:  Tg  =  3000,  5000,  and  8000°H 

(e)  Carrying  out  plume-induced  separated  flow  calculations.  The  resulting  con¬ 
servation  equations  for  total  mass  and  species  fluxes  for  different  Xg  are 
illustrated  in  Figures  1-5  through  1-7.  The  intersection  points  satisfy  both 
mass  and  species  conservation  requirements. 

(f)  Cross-plotting  Tg  versus  Yg  (indicated  by  BLSEP)  and  intersecting  with 
the  corresponding  curve  based  on  equilibrium  thermochemical  calculations 
(NASA-Lewis).  The  intersection  points  now  satisfy  mass,  species,  and 
energy  conservation  requirements  (Fig.  1-8).  Thus,  the  proper  Yg  values 
(fuel-air  mixture  ratios  in  the  base  region)  are  determined. 
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Figure  1-6  Conservation  of  Mass  and  Species  Fluxes  for  Xj 
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(g)  The  corresponding  discriminating  streamline  locations,  Vs  are  ascertained 
with  additional  cross-plot  (Fig.  1-8). 


(h)  Next,  the  discriminating  streamline  total  pressures  are  plotted  for  various 
values  of  nse  and  Yp  (Fig.  1-9),  resulting  in  the  determination  of  the  pro- 
per  Pf  /Pn  values  for  different  xs. 
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(i)  Independent  of  the  separated  flow  calculations,  the  trailing  shock  pressure, 

Pw,  formed  by  the  intersection  between  the  stream  sin  ar  layer  and  the  plume 

shear  layer  can  be  determined  on  the  basin  of  the  classic  wave  interaction 

concept  (Sec.  n.  Par.  5. a).  The  interpolation  scheme  and  the  resulting  P  for 

w 

xg  =  20,  30,  and  40  feet  are  depicted  in  Figure  1-10.  Using  a  Nash  factor  of  0.56, 
the  following  results  are  obtained: 


Calculated 


^w 

Els. 

xs 

(ft) 

PB 

PB 

20 

9.6 

5.8 

30 

10.2 

6.15 

40 

10.9 

6.55 

(j)  Now  cross -plotting  the  Ptg/Pg  values  base'i  on  the  shock  interaction  theory, 
a  final  Pts/pg  is  obtained.  Inasmuch  as  tl  is  Ptg/Pp  satisfies  all  the  con¬ 
servation  requirements  as  well  as  the  reatta*  hment  criterion,  it  represents 
the  solution.  The  corresponding  Xg  is  found  to  be  38.5  ft  from  the  tip. 

The  plane -of- symmetry  solution  delineated  above  represents  the  basic  plume-induced 
separated  flow  solution.  Pseudo-two-dimensional  solutions  for  other  strips  and  a 
special  technique  to  construct  the  separation  footprint  have  been  discussed  in  detail  in 
par.  1  of  Section  El  and  will  not  be  repeated  here. 
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Figure  1-9  Streamline  Total  Pressures  for  Various  VSq  and  YR 
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Figure  I- 10  Determination  of  Trailing  Shock  Pressure  and  Separation  Distance 
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Appendix  II 

LAMINAR  CROSS-FLOW  BOUNDARY  LAYER  ANALYSIS 

by  M.  C.  Fong 

The  flow  Held  surrounding  a  maneuvering  lifting  reentry  vehicle  is  three-dimensional 
not  only  due  to  the  body  geometry  but  also  due  to  the  attitude  and  orientation  of  the 
vehicle  under  hypersonic  flight  conditions.  Peculiar  surface  streamline  patterns  have 
been  observed  (Refs.  37,  53,  and  54)  but  analytical  prediction  of  such  streamline  devia¬ 
tion  behavior  and  the  concomitant  viscous  flow  properties  has  not  reached  beyond  a 
semiempirical  stage  except  for  isolated  Instances  (Ref.  55).  However,  the  local  in- 
viscid  and  boundary  layer  flow  properties  must  be  well  understood  in  order  to  obtain 
a  satisfactory  plume-induced  separated  flow  solution  presented  in  this  analysis.  De¬ 
velopment  of  an  analytical  method  to  predict  the  attached  flow  properties  under  realis¬ 
tic  LRV  flight  conditions  thus  constitutes  the  main  goal  of  this  appendix. 

For  expediency,  a  strip  method  has  been  adopted  for  both  attached  and  separated 
boundary  layer  predictions  in  conjunction  with  the  plume-induced  separation  problem. 
The  strip  method  essentially  rationalizes  that  the  two-dimensional  flat-plate  solution 
within  each  narrow  chordwise  strip  is  approximately  valid  and  any  Influence  from  the 
adjacent  strip  is  of  secondary  Importance.  However,  the  streamline  deviation  cited 
above,  which  is  generally  attributable  to  inviscid-viscous  interaction  at  hypersonic 
speeds,  necessitates  a  re-examination  of  the  strip  method  in  application  to  the  LRV 
maneuvering,  or  at  least  establishment  of  a  set  of  criteria  within  which  the  strip 
method  is  satisfactory  for  predicting  the  exhaust  plume  effects . 

In  this  analysis,  an  engineering  method  is  proposed  to  resolve  the  streamline  deviation 
problem,  in  that  the  cross-flow  boundary  flow  in  terms  of  the  streamline  coordinate 
system  is  treated  with  the  small  perturbation  approach.  Then,  using  an  available  sur- 

I 

face  oil  flow  pattern,  the  unknown  lnviscid  streamline  is  determined  which  can  be 
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compared  with  the  corresponding  streamline  based  on  the  strip  method.  This  inverse 
approach  avoids  the  displacement  thickness  growth  problem  and  the  subsequent  inviscid- 
viscous  flow  interaction  phenomenon  and  therefore,  permits  an  assessment  of  the  strip 
method.  The  predicted  inviscid  streamline,  of  course,  provides  an  input  for  plume- 
induced  separated  flow  prediction. 

1.  FORMULATION  OF  GOVERNING  DIFFERENTIAL  EQUATIONS. 


The  formidable  task  involved  in  the  three-dimensional  boundary  layer  problem  lies  in 
the  cross-flow  generation  and  streamline  distortion  within  the  boundary  layer  for 
which  the  high-order  nonlinearity  generally  precludes  any  similarity  treatment  or 
other  simplifications.  However,  by  using  the  streamline  coordinate  system  for  which 
one  of  the  coordinates  is  set  to  be  coincident  with  the  local  outer-edge  streamline  pro¬ 
jected  in  the  tangent  plane,  it  can  be  envisioned  that  the  governing  equations  would  be 
subject  to  small  perturbation  treatment  except  for  large  curvatures  in  the  flow  field. 
This  analysis  essentially  follows  the  formulation  presented  in  References  41  and  48, 
and  places  emphasis  on  development  of  a  small  cross-flow  solution  that  can  be  easily 
incorporated  in  the  plume-induced  separated  flow  analysis. 


Similar  to  the  consideration  by  Beckwith  (Ref.  47),  a  three-dimensional  orthogonal 
streamline  coordinate  system  (x,y,z)  is  introduced  with  corresponding  velocity 
component  u  tangent  to  the  external  streamline,  v  normal  to  the  body  surface,  and 
w  the  cross-flow  velocity  component  normal  to  u  in  the  tangent  plane  (Fig.  II— 1) .  * 
The  length  elements  are  defined  as: 


ds  =  (x  ,  z)  dx 
dy  =  dy 

dn  =  e9  (x  ,  z)  dz 


(0-1) 


*The  coordinate  system  is  left-handed;  the  equations  are  identical  to  those  by  Beckwith 
(Ref.  47)  if  the  coordinates  y  and  z  and  the  velocity  components  v  and  w  are 
interchanged. 
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Figure  II- 1  Streamline  Coordinate  System 


and  the  outer-edge  conditions  require  that 


u  =  ue 

w  =  w  =  0 
e 


(D-2) 


With  the  assumptions  of  perfect  gas  and  the  absence  of  diffusion  and  chemical  reactions, 
the  governing  equations  for  a  three-dimensional  laminar  compressible  boundary  layer 
are  then  expressible  as 


♦ 
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Continuity 


Tj  I  <pue2»  +  If  (pvc2»  *  T[  H  <pvei>  '  0 


(II-3) 


x-Momentum 


Be  2  Be 

u  3u  A  „  3u  .  pw  8u  ,  puw  1  p w  2  1  dp 

pv-5y  e^5z  IT  '  W  = 


(0-4) 


y-Momentum 


If  '» 


(0-5) 


z -Momentum 


Oa  O  Qn 

u  1w  x  .  8w  pw  8w  puw  2  pu  1 

pe'1  Jx  pv3y  ~e^  "Sz" +  e,e„  8x  ~  e,e„  8z 


'12 


12 


■  --s:  If +  w  ("if)  <n'6) 


Energy 


0  u  9H  .  flv  9H  .  pw  3H  =  8  [  /8H  1  -  Pr  9h\|  m__. 

p  ej  8x  pv  8y  e2  8z  8y  r\57  Pr  8y/l  * 


From  the  continuity  equation,  two  stream  functions  ip  and  p  can  be  defined  as  follows: 


pue 


8  ip  ( dip  30\  1  30 

2  =  87  *  pVe2  =  -fe+8F/i;  ’  pwei  =37 


(0-8) 


Although  for  small  perturbation  approach,  these  partial  differential  equations  can  be 
simplified  by  examining  the  crossflow  terms  (Refs.  47  and  48),  the  present  analysis 


159 


follows  the  scheme  discussed  in  Reference  56  so  ‘hat  a  system  of  similarity  equations 
in  successive  order  can  be  obtained.  To  accomplish  this,  the  following  transforma¬ 
tion  variables  are  introduced: 


x  s 

4  =  /  p  u  e  dx  =  /  puds 

J  Wo  1  /  Wo 

o  o 


(II-9) 


where  subscript  o  denotes  some  reference  value,  such  as  that  for  local  stagnation 
conditions.  Furthermore,  the  stream  functions  are  expressed  in  the  following  manner: 


It  follows  that 


♦  =  /2Tue  e2  f(£>  £) 


0  =  y/^e  eiS(l>n>£) 


(II- -10) 


u  =  u  and  w  =  u  -§§■ 

e  arj  e  dv 


(11-11) 


Also,  at  the  outer  edge  of  the  boundary  layer 


„  u  9u 
pe  e  e 

e,  ax" 


1  9p 

ej  ax 


(II- 12) 
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and 


du 

e 

az 


o  u 

He  e 


3P 

9z 


(n-13) 


Substituting  these  transformation  variables  into  the  governing  Equations  (II-3)  to  (II-7) 
and  after  algebraic  rearrangement,  the  following  equations  are  obtained: 

Streamwise  Momentum 


a  Lit} 
®i  X  a 


~2 

an 


g  a  f  af  ag 

?  dr}2  dv  3tj 
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Cross-Flow  Momentum 
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(U-14) 
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And  with  the  definitions  of  the  total  enthalpy 


e 


(11-16) 


and 


(11-17) 


the  energy  equation  becomes: 
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♦  M  !fi  tgH) 
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=  0 


(11-18) 


where  N  =  Pp/p0M0  is  a  viscosity  function  which  can  be  simplified  to  become  a 
pressure  ratio  N  =  P/PQ  if  a  linear  viscosity-temperature  relation  holds. 
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The  boundary  conditions  are: 


At  the  wall,  tj  =  0 

u  =  w  =  0  :  =  0  ,|i=0 

dV  dV 

v  =  0  :  f  =  0,  g  =  0 
At  the  outer  edge,  rj  -•  <*, 


u  =  u  ,w  =  w  =  0,H  =  H 
e  e  e 
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(11-19) 


Introducing  the  following  definitions: 
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Streamline  pressure 
gradient  parameter 


Curvature  parameter 


Dilatation  parameter 


Cross-flow  pressure 
gradient 
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(11-20) 


Equations  (U-15),  (U-16),  and  (11-18)  become  (Ref.  56) 
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2.  PERTURBATION  GENERATION. 


Small  cross  flow  is  attained  when  the  surface  curvature  and  dilatation  parameters  as 
well  as  the  cross-flow  pressure  gradient  are  small  in  comparison  with  the  stream- 
wise  velocity  component.  Under  this  circumstance,  a  small  parameter,  e  ,  may  be 
introduced  such  that 
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(11-24) 


where  the  subscripts  0,1  , 2  . . .  represent  the  zero-,  first-,  second-order  solutions, 
etc.  For  simplicity,  the  following  assumptions  are  also  introduced: 

•  All  the  f,  g,  and  0  functions  are  similar,  i.e. ,  dependent  on  t?  only. 

•  Prandtl  number  equal  to  unity 

•  The  viscosity  parameter  N  can  be  approximated  by  unity,  this  being  admissi¬ 
ble  locally  for  relatively  flat  surfaces  considered  herein. 

Equations  (11-21)  to  (11-23)  then  become: 

Zero-Order 


fM  +  {  { ...  A  +  0_\  +  p/g  _  f.2\  =  0 

o  o  o  \  2m '  \  0  o  / 


e"  +  f  0'  (\  +  -4-)  =  o 

o  oo\  2a/ 


(11-25) 

(11-26) 
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First-Order 
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And  so  on.  Note  that  (3 *  -  -K^  is  utilized  in  obtaining  Equation  (11-28) . 

For  hypersonic  flow,  these  equations  can  be  further  simplified  because  (3/i 
negligibly  small.  Then,  we  obtain 

Zero-Order 


f'"  +  f  f"'  +  B  (e  -  f'2) 

o  o  o  Ho  o  / 


=  0 


0"  +  f  0' 
n  0  0 


=  0 


First-Order 


fr +  vr  -  <2w +  fofi +  +  f0foK2i> 


=  0 


g;M  -  w  +  h  i9o  - f;2 )  =  0 


ev  +  f  o;  +  eif,  +  f  9'K„  =  o 

1  01  01  0021 


where  (3*  =  /?*m  is  a  zero-order  quantity  as  seen  from  the  definition  of  (31 


(11-27) 

(11-28) 

(11-29) 

becomes 

(11-30) 

(H-31) 

(11-32) 

(11-33) 

(11-34) 
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It  is  seen  that  equations  (11-30)  and  (11-31)  are  identical  to  the  transformed  equations 
obtained  by  Cohen  and  Reshotko  (Ref.  57)  if  6  is  replaced  by  S  with 

S  =  Q  -  1  (11-35) 

Furthermore,  the  first-order  cross -flow  Equation  (U-33)  depends  on  f  and  Qq  but 
is  independent  of  the  other  first-order  functions,  f^  and  0^  .  Hence,  if  attention  is 
focused  on  the  solution  of  Equation  (11-33),  fj  and  flj  need  not  be  solved  at  all.  It 
should  be  noted  also  that  because  the  f  and  dQ  values  are  already  available,  Equa¬ 
tion  (11-33)  is  linear  in  gj  for  which  the  numerical  solution  is  quite  easy  to  obtain 
despite  a  two-point  boundary  condition  involved. 

4.  COMPUTATION  SCHEME. 

The  widely  used  Runge-Kutta-Gill  forward  integration  technique  was  employed  to 

solve  Equa*‘?n  (11-33)  and  a  computer  program  for  UNIVAC  1108  was  developed  to 

carry  out  the  numerical  calculations.  The  zero-order  values  f  ,  9  ,  f'  etc.  were 

o  o  o 

obtained  from  the  Cohen-Reshotko  results  (Ref.  57)  which  were  then  used  as  part  of 
the  input  data.  The  problem  concerning  the  two-point  boundary  condition  was  re¬ 
solved  by  replacing  the  outer-edge  condition  by  an  assumed  wall  condition  with  which 
a  tentative  cross-flow  solution  was  obtained.  Following  comparison  with  the  pre¬ 
scribed  outer-edge  condition,  a  new  wall  condition  was  then  assigned.  Iteration  on 
this  particular  boundary  condition  would  eventually  lead  to  the  desired  solution. 

5.  RESULTS  AND  DISCUSSION. 

A  number  of  cases  in  terms  of  different  0  and  Sw  (wall  temperature  parameter) 
have  been  investigated.  The  results  are  presented  in  Figures  n-2  through  11-11  which 
show  the  cross-flow  velocity  profiles  for  different  values  of  0*  .  In  several  cases, 
the  result  exhibits  a  reverse  profile,  but  with  only  a  few  cases  at  hand,  parameters 
governing  the  reverse  flow  could  not  be  readily  established. 
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The  reason  that  0*  values  were  selected  to  be  of  the  order  of  10  is  based  on  the 
flow  conditions  used  in  Reference  49,  i.e.,  M  =  6.8,  Tt  =  150(TR  .  Fora 
6 -degree  oil  flow  streak  shown  in  Figure  45,  is  estimated  as  follows: 
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For  a  flow  deviation  of  6  degrees,  8ue/8z  is  found  to  be  0. 105.  It  is  assumed  that 
this  value  would  be  that  for  e  also,  from  which  an  approximate  value  of  0.02  for 
/?*  results.  As  a  comparison,  the  Blasius  profile  which  typifies  the  main  stream 
velocity  distribution  is  about  an  order-of-magnitude  greater  than  the  cross  -flow 
profile  (Fig.  11-12).  Hence,  under  the  small  perturbation  conditions,  the  cross-flow 
effect  is  indeed  of  secondary  importance. 


The  cross-flow  velocity  profiles  can  be  utilized  for  wall  shear  and  heat  transfer  esti¬ 
mates.  However,  the  main  concern  here  lies  in  the  determination  of  inviscid  stream¬ 
line  patterns,  inversely,  based  on  the  cross-flow  solution;  the  other  aerothermodynamic 
aspects  mentioned  above  have  not  been  explored. 


Results  of  inviscid  streamline  prediction  and  comparison  with  the  strip  theory  for  the 
sample  delta  surface  are  discussed  in  Section  II  and  depicted  in  Figure  45. 
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Figure  II -3  Cross-Flow  Velocity  Profile  for  (3  =  -0.1  and  S  =  -0.8 
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Figure  n-8  Cross-Flow  Velocity  Profile  for  0  =  0.5 
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Figure  II-9  Cross-Flow  Velocity  Profile  for  0=  1.0  and 
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Figure  11-11  Cross-Flow  Velocity  Profile  for  /?  =  2  9  and 


APPENDIX  III 


This  section  presents  illustrations  for  predicted  separation 
footprints  on  both  the  windside  and  leeside  surfaces  of  the 
simple  lifting  body  configuration  considered  in  the  text.  A 
total  of  43  illustrations  with  various  Mach  numbers,  altitudes, 
and  angles  of  attack  is  included  in  this  Appendix;  however, 
those  which  display  negligibly  small  separation  footprints  are 
excluded. 

The  main  physical  parameters  are: 

(  =  chord  length  =  45  ft 

=  slope  length  =  f  sec  10° 

b  =  Windside  base  length  *  15.6  ft 

b1  =  leeside  base  length  =  0.707  b 

d  =  Centroid  location  measured  from  base 
c 

A  =  Separation  footprint  area 
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A  »  .053  b/ 

d  *  .093/ 
o 


Figure  m-5  Separation  Footprint  on  Windside  Surface  for 
■  6,  280-Kft  Alt,  and  a  «  0  deg 
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Figure  HI-22  Separation  F< 
M03=  10,  23 
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Figure  m-36  Separation  Footprint  on  Leeside  Surface  for 
M*  =  10,  230-Kft  Alt,  and  a  =  +20  deg 
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A  modified  version  of  the  Chapman-Korst  model  involving  mass  diffusion, 
equilibrium  thermochemistry,  and  appropriate  semiempirical  relations  is 
utilized  to  predict  the  overall  flow  characteristics  within  the  plume-induced 
separated  flow  region  during  synergetic  maneuvering  of  a  lifting  reentry  ve¬ 
hicle.  The  exhaust  plume  selected  in  conjunction  with  the  separated  flow 
analysis  is  predicted  on  the  basis  of  the  method-of-charaeteristics  solution 
with  an  average  spccific-heat  ratio  of  1.38.  Thermochemical  equilibrium 
is  assumed  to  exist  within  the  plume-induced  separated  flow  region,  and 
equilibrium  temperature  curves  for  various  altitudes,  vehicle  Mach  num¬ 
bers,  and  propellant-air  mixture  ratios  are  presented.  The  strip  method 
is  incorporated  in  the  analysis  to  estimate  the  three-dimensional  separated 
flow  characteristics  caused  by  the  vehicle  geometry  and  vehicle  attitude 
during  the  maneuver.  The  extent  of  validity  of  the  strip  method  is  estab¬ 
lished  with  the  aid  of  a  cross-flow  boundary  layer  solution.  On  the  basis 
of  the  present  analytical  model,  the  propulsion  effects  on  aerodynamic 
characteristics  of  a  simple  lifting  body  configuration  with  delta  planform 
and  triangular  cross  section  under  a  wide  spectrum  of  freestream  condi¬ 
tions  have  been  evaluated  and  illustrated  by  appropriate  design  charts. 
Although  the  usefulness  of  the  prediction  method  has  been  demonstrated, 
additional  relevant  data  are  required  to  finalize  the  theory  development 
and  to  ensure  a  broad  applicability  of  the  analytical  model. 
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ERRATA 


In  the  AFFDL  technical  report  entitled  "Propulsion  Effects 
on  Aerodynamic  Characteristics  of  Lifting  Reentry  Vehicles", 
by  Michael  C.  Fong  and  Carl  F.  Ehrlich,  Jr.  (AFFDL-TR-70-12 , 
March  1970),  the  following  errors  should  be  corrected. 


Page  xvi,  the  7th  line  from  the  bottoms  The  definition  of  A 
should  be  /8*  ^ 

/  1° 

Page  135 ,  line  8s  — —  should  be  changed  to 

Page  138,  equation  (I -1*0):  The  integrand  should  be  multiplied  by  jP. 

Page  139,  equation  (1-1*6):  The  integrand  should  be  multipliea  by^2  . 

Page  11*0,  equation  (1-1*7):  The  integrand  should  be  multiplied  by 

However,  13  c  was  introduced  in  an  alternate  form  of 

energy  equation  and  is  not  relevant  to  any  of  the  equations 
used  in  this  report. 

Page  U*0,  equation  (I-U9):  The  lower  limit  of  the  integral  should  be 
y  instead  of  y  ,  and  the  "negative"  sign  in  functions 

Se  _  3j 


Page  11*1,  equation 


F^(-  tsQ  1  CP  )  and  F^  C -  Vse>°°  )  should  be  removed. 

equation  (1-50):  The  "negative"  sign  in  functions  ,<©) 

and  >  00  )  should  be  removed.  sj 


Page  li*l,  equation  (1-51):  The  equation  should  be  corrected  to 

X  +  Y  =  =  f^udy 

J  \  J*j  JAj 

Page  11*1,  equation(I -52 )  s  The  equation  should  be  corrected  to 

y~ y  -  JZs1  v  f  r/’/  oo'i  -  r  / /  _  *Quiyi  T m 


j  j 

+  ^  [%!,%•,»)  +  0  -r8)i2^j»;J  =  o 

J 

Page  llx2f  equation  (I-55)i  The  equation  should  be  corrected  to 

jy  =  r,  [ 

+  [yb  1,(  %.»>  +( 1  -Y.) 

J 


